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Abstract 

The origin of the hot phase of the early universe remains so far an unsolved puzzle. 
A viable option is entropy production through the decays of heavy Majorana neutri- 
nos whose lifetimes determine the initial temperature. We show that baryogenesis 
and the production of dark matter are natural by-products of this mechanism. As 
is well known, the cosmological baryon asymmetry can be accounted for by lepto- 
genesis for characteristic neutrino mass parameters. We find that thermal gravitino 
production then automatically yields the observed amount of dark matter, for the 
gravitino as the lightest superparticle and typical gluino masses. As an example, 
we consider the production of heavy Majorana neutrinos in the course of tachyonic 
preheating associated with spontaneous B — L breaking. A quantitative analysis 
leads to constraints on the superparticle masses in terms of neutrino masses: For 
a light neutrino mass of 10~ 5 eV the gravitino mass can be as small as 200 MeV, 
whereas a lower neutrino mass bound of 0.01 eV implies a lower bound of 9 GeV on 
the gravitino mass. The measurement of a light neutrino mass of 0.1 eV would rule 
out heavy neutrino decays as the origin of entropy, visible and dark matter. 
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1 Introduction 



Detailed studies of the cosmic microwave background provide direct evidence for the hot 
thermal universe close to its minimal temperature pQ. The extrapolation to higher tem- 
peratures, beyond primordial nucleosynthesis, is very uncertain, and we do not know 
how large the maximal temperature of the hot early universe has been. It is widely be- 
lieved that the universe was 'reheated' by a transition from a preceding inflationary phase 
where 'vacuum energy' dominated the expansion [2J. Knowing the resulting reheating 
temperature is of fundamental importance since it is closely related to the origin of the 
matter-antimatter asymmetry and the nature of dark matter. 

In a recent paper we have suggested that the entropy of the hot early universe was 
produced in the decays of heavy Majorana neutrinos whose lifetimes determine the re- 
heating temperature (3j. We have demonstrated that the baryon asymmetry and dark 
matter then naturally result as by-products: The mechanism of baryogenesis is a mixture 
of thermal and nonthermal leptogenesis, and the dominant component of dark matter is 
the gravitino which is assumed to be the lightest superparticle. 

Our work is closely related to previous studies of thermal leptogenesis [4"f5] and non- 
thermal leptogenesis via inflaton decay [6117] , where the inflaton lifetime determines the 
reheating temperature. In supersymmetric models with global B — L symmetry the scalar 
superpartner N± of the lightest heavy neutrino N\ can play the role of the inflaton in 
chaotic [8] or hybrid [9] inflationary models. Coherent N\ oscillations after inflation may 
even dominate the energy density of the universe [10]. Nonthermal leptogenesis can be 
realized at significantly smaller temperatures than thermal leptogenesis. In this way the 
'gravitino problem' for heavy unstable gravitinos [TTfl2|ll3| can be easily avoided, which 
has been one of the main motivations of nonthermal leptogenesis. 

It is well known that the high temperatures characteristic for thermal leptogenesis 
can become a virtue if the gravitino is the lightest superparticle (LSP). For superparticle 
masses as they arise in gravity or gaugino mediation, thermal production of gravitinos can 
then explain the observed amount of dark matter [14]. As pointed out in [3], the required 
high temperatures are indeed realized if the universe is reheated through the decays of 
the heavy Majorana neutrinos. The vacuum decay width of the lightest heavy Majorana 
neutrino is given by 



where we have used the typical values Mi ~ 10 GeV for the N\ neutrino mass, 
rh\ ~ 0.01 eV for the effective light neutrino mass and t> EW = 174 GeV for the vac- 




(1) 
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uum expectation value of electroweak symmetry breaking. The corresponding reheating 
temperature is given by 

90 ^ 1/4 



T ™~{8^Z) ^%M P ~lV°GeV, (2) 

where we have used g+ tP ~ 200 for the effective number of relativistic degrees of freedom, 
and M P = 1.22 x 10 19 GeV is the Planck mass. 

From Eqs. and (|2| one obtains the reheating temperature Trh in terms of the 
neutrino masses rhi and Mi, the two key parameters for thermal and nonthermal lepto- 
genesis. Assuming the gluino to be the heaviest gaugino, thermal gravitino production 
is dominated by QCD processes, which yields the gravitino abundance in terms of the 
gravitino mass mg and the gluino mass rrig, 

where the coefficient C = 0.26 to leading order in the gauge coupling [T5lfT6] |^] Since Trh 
depends on mi and Mi, the requirement^] fig h 2 = Odm'* 2 — 0.11 [20j yields a connection 
between neutrino and superparticle mass parameters. The neutrino masses rhi and M\ 
are in turn constrained by the condition that the maximal baryon asymmetry is larger 
than the observed one, 7]b > r]° B hs = 6.2 x 10~ 10 [20]. 

In the following sections we shall study in detail the connection between neutrino 
and superparticle masses, which is implied by successful leptogenesis and gravitino dark 
matter. As an example, we shall consider tachyonic preheating [21], associated with 
B — L breaking, as a mechanism which can lead to a phase where the energy density is 
dominated by heavy Majorana neutrinos. As we shall see, the final baryon asymmetry 
and the dark matter abundance can then be calculated in terms of several parameters of 
the Lagrangian, independent of initial conditions: the scale vb-l of B — L breaking, the 
heavy Majorana neutrino mass M 1; the effective light neutrino mass rhi, the gravitino 
mass mg and the gluino mass m~ g . Particularly interesting is the resulting connection 
between the lightest neutrino mass mi and the gravitino mass. 

Our analysis requires a flavour model which is flexible enough to allow for a large range 
of the neutrino masses Mi and rhi, the crucial parameters for leptogenesis. Such a model 
is described in Section 2. Subsequently, in Section 3, we discuss in detail the Boltzmann 

^ote that C has anO(l) uncertainty due to unknown higher order contributions and nonperturbative 
effects [TS]. Resummation of thermal masses increases C by about a factor of two [T7] . 

2 For the superparticle masses considered in this paper the contribution to the gravitino abundance 
from the decay of the next-to- lightest superparticle (NLSP) is negligible |18|19j . 



equations that describe the time evolution of the state produced in the tachyonic decay of 
the false vacuum with unbroken B — L symmetry. A novel technical aspect of our analysis 
is the separate treatment of thermal and nonthermal contributions to the abundance of 
the heavy Majorana neutrinos. 

Using the set of Boltzmann equations derived in Section 3, we study an illustrative 
example of our mechanism in Section 4. The parameters are chosen such that the non- 
thermal contribution to leptogenesis dominates. Particular emphasis is given to the time 
dependence of the various production and decay rates and the emergence of a plateau 
where the temperature is approximately constant. In Section 5 the analysis is extended 
to the entire parameter space. First, the ranges of M\ and fhi are determined for which 
leptogenesis is successful. The observed dark matter abundance then constrains the su- 
perparticle masses in terms of the neutrino masses. 

Our results are summarized in Section 6. The appendices deal with various technical 
aspects of our calculations: conventions for the Boltzmann equations in Appendix A, the 
distribution function of thermally produced neutrinos in Appendix B, analytical approxi- 
mations for the reheating temperature in Appendix C and semi-analytical results for the 
gravitino abundance in Appendix D. 

2 Flavour model and leptogenesis 

In the following we shall describe a flavour model which describes masses and mixings of 
quarks, charged leptons and neutrinos, and which is flexible enough to allow for a large 
range of the neutrino parameters Mi and rh\ that are crucial for leptogenesis. The model 
is a variant of [22] and satisfies all constraints from flavour changing processes [23J . 

We consider the extension SU(3) c xSU(2)ixU(l)y xU(1)b_i of the standard model 
gauge group. The Yukawa interactions of quarks and leptons with Higgs fields are de- 
scribed by the following superpotential for matter superfields, 

W M = /^lOilOjiZ,, + hp*10 3 H d + !>;/>] n'M,, + ^n^S . (4) 

Here the standard model fermions have been arranged in SU(5) multiplets, 10 = (q, u c , e c ) 
and 5* = (d c , £), and i,j — 1 . . .3 are flavour indices. n c contain u R , the charge conjugates 
of the right-handed neutrinos, which are related to the heavy Majorana neutrinos iV 
through iVj = v Ri + v c Ri . For simplicity, we have used SU(5) notation assuming that the 
colour triplet partners of the Higgs doublets have been projected out. 

Vacuum expectation values of the scalar Higgs fields, (H U) d) = v Ut d, break the elec- 
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Table 1: Chiral U(l) charges. 



troweak symmetry. The superpotential 

W B „ L = (v%_ L - 2S 1 S 2 ) (5) 

enforces B — L breaking, with (Si) = (S 2 ) = vb-l/V%, via the supersymmetric Higgs 
mechanism. Using the field b asis Si — S exp (i$) / y/2, S 2 = S' exp (— i$) / ^/2 and shifting 
around the vacuum expectation value, 5" = Vb-l + S, S and T have a common Dirac 
mass term, whereas $ and the B — L vector multiplet form together a massive vector 
multiplet. Since B — L is gauged the inflaton is identified with the scalar component of 
the singlet field T. The potential for the scalar neutrinos is not sufficiently flat. 

The pattern of Yukawa couplings is determined by a Froggatt-Nielsen U(l) flavour 
symmetry, following Ref [22]. The matter fields tpi and an extra singlet S carry charges 
Qi and —1 under the flavour symmetry respectively, and are coupled together via effective 
non-renormalisable interactions associated with a scale A > Aqut- The Yukawa couplings 
are generated once the flavour symmetry is spontaneously broken by the expectation value 
of the S field, and are given by 

hij oc r] Q > +Q J , (6) 

where 77 = (X) / A. The hierarchies of quark and lepton masses are then naturally obtained 
for rj 2 ~ 1/300 using the chiral charges listed in Table [TJ It is important to note that the 
Yukawa couplings are only specified up to factors of 0(1). Cosmology further constrains 
the chiral charges. For instance, successful thermal leptogenesis requires a + d = 2 [22J. 

We shall restrict our analysis to the case of a hierarchical heavy neutrino mass spec- 
trum, Mi M 2) 3, which is obtained for flavour charges of Ni and ^2,3 separated by 
one unit, b = c = d — 1. This is sufficient to illustrate our main point, the contraints 
imposed on the gravitino mass by neutrino masses. The lepton asymmetry will mostly 
be generated by decays of the lightest heavy neutrino N%. The masses of the Majorana 
neutrinos are given by 

M x ~ v 2d v B _ L , (7a) 
M 2 , 3 ^ V 2(d ' 1] vb-l , (7b) 



6 



where M 2 /M 3 = 0(1). The mass spectrum is now parametrised by the three remaining 
free charges a, d and e, which can be related to the physical parameters vb-l, Mi and 
the Higgs mass ms- 

Using the Dirac mass matrix mo = h u v u and the Majorana mass matrix M = h n VB-L 
derived from Eq. (|4]), the eigenvalues of the light neutrino mass matrix m = —mvM~ 1 m T D 
are given by 

m x ~ ^«+2^|w (8a) 

V B -L 

m 2 ~m 3 ~ r] 2a , (8b) 

v B -L 



where we have used tan/3 = v u /vd > 0(1) and v u ~ t>EW- Choosing m v = ^Jm^H — 
\(m\ — m^)(m\ — m 2 )| 1//4 ~ 3 x eV (cf. PQ), the B — L breaking scale is determined 
by the flavour charge a, 



VB-L^Tp^-, (9) 

where u| w /m„ ~ 10 15 GeV is the grand unification mass scale. Once vb—l is fixed, Mi 



is directly related to the charge d through Eq. (7a). 

The ranges over which the chiral charges a and d, and thus the physical parameters 
v b-l and Mi, are allowed to vary is restricted. First, the requirement that Yukawa 
couplings do not exceed the top- Yukawa coupling imposes the lower bounds a > and 
d > 1. Furthermore, the upper bound a < 1 follows from tan ft > 0(1). No corresponding 
upper bound on the charge d exists but, as we shall see later, a CP asymmetry sufficiently 
large for successful leptogenesis requires d < 3. Using Eq. ([9]), the allowed range of B — L 
breaking scales reads 

3 x 10 12 GeV < v B ~L < 1 x 10 15 GeV . (10) 



For fixed B — L breaking scale, the possible range of Mi is given by Eq. (7a) 



1 x 10 5 GeV < Mi < 1 x 10 10 GeV for v B -l = 3 x 10 12 GeV , (11a) 
3 x 10 7 GeV < Mi < 3 x 10 12 GeV for v B -l = 1 x 10 15 GeV . (lib) 

The ranges for vb-l an d Mi correspond to a continuous variation of the flavour charges, 
which can effectively be realised by fractional charges. Note that the constraint from 
thermal leptogenesis, a + d = 2, is now relaxed. 
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Given the Yukawa couplings, one easily obtains the decay widths of the heavy Majo- 
rana neutrinos, 

n rhj Mf , , 

* . (12) 



with the effective light neutrino masses 



1 v 2 
m i = T7 ( m D m D)u - V 2a - m u • (13) 

M Vb-L 



The CP asymmetry in the heavy neutrino decays are given by [241125] 

e i= - 1 — yimiUh^h").] 2 ] f(^A , (14) 

BK(h^h v ) u ^ U v h n J \MiJ K J 

where we use the standard model expression for Using our flavour model, one gets 

e x ~ 0.1 t/ 2 ^ = 0.1 ^1 , e 2i3 ~ Cl ^ 2 • (15) 

Note that this is the maximal CP asymmetry for fixed Mi [TUll2"6"] , which is obtained in the 
limit fhi -> 0. For Mi ~ 10 10 GeV, this yields ei ~ 10 6 . For other Majorana neutrino 
masses, the asymmetry scales like the mass ratio, e% ~ 10 _6 Mi/10 10 GeV. 

Since the light neutrino mass matrix is not hierarchical, the 0(1) uncertainties in the 



h p Yukawa couplings can lead to large deviations from the relation ( 13 ) between fhi and 
m v . The only rigorous inequality is fhi > mi [27]. We take these uncertainties into 
account by varying the effective neutrino mass in the range 

10~ 5 eV < rhi < 0.1 eV . (16) 

Since the heavier Majorana neutrinos only play a marginal role in our scenario (see 



below), we ignore possible deviations from the relation (13) and use 777-2,3 = "V- 

In the following we consider the 'waterfall transition' from the false vacuum (S) = 
to the true vacuum (S) = vb~l, which may happen at the end of hybrid inflation. A 
tachyonic instability in the Higgs potential leads to spinodal growth of the long- wavelength 
Higgs modes. The true vacuum is reached after a rapid transition at time tpn [28J, 

<^)U„ = ipH ,J- ln «. (17) 
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The expression in the supersymmetric standard model would only slightly increase the value of ej. 
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The Higgs boson forms a massive supermultiplet together with a Dirac fermion and three 
additional bosons. Its mass and the false vacuum energy density are given by 

m 2 s = \v\_ L , po = \ Xv b-l • ( 18 ) 

The size of the coupling A is determined by the flavour charge e in Table 1. We shall 
restrict our analysis to the case e = 2(d — 1), such that mg — M2.3. As a consequence, the 
Higgs boson only decays to pairs of Ni neutrinos and not to pairs of N 2 or N 3 neutrinos. 
Like the heavy neutrino masses and the CP asymmetries, the Higgs mass only depends 
on the flavour charge a + d whereas the false vacuum energy density is determined by a 
different combination of charges, 



ms „ 2(o+d-l) 1/4 „ 2a+d-l ^V (lg) 

m„ ' m„ 



During the tachyonic preheating the energy of the false vacuum is converted mostly 
into a nonrelativistic gas of S bosons {\ps\/ m s 1), with an admixture of heavy neu- 
trinos. Their contribution to energy density and number densities is determined by their 
coupling to the Higgs field [28J, 

r W , = -^1.5x 10- 3 g N Xf(a h 0.8) , (20a) 
Po 

n Ni ~ 3.6 x 10" 4 g N m% f{a, 0.8) /a , (20b) 

where = 2 and f(a, 7) = \J a 2 + 7 2 — 7 with on = hf/VX- For the heaviest neutrinos 
A^2,3, one obtains 

r^-lxlO- 3 ^) 2 , (21) 

while the relative contribution from the lightest right-handed neutrinos, r^/r^^ ~ 
0(r] 2 ), is negligible. 

For simplicity, we neglect all superpartners as well as B — L gauge bosons and inflaton 
modes (T) which are produced during tachyonic preheating. We expect their contributions 
to not significantly change our results, similarly as in supersymmetric leptogenesis [29] . 
A detailed discussion will be presented in [30J . A further important aspect of tachyonic 
preheating is the production of cosmic strings [21] ■ Their effect on the baryon asymmetry 
is model dependent |31j. It has to be analysed for the parameters of our model [30] taking 
also into account non-minimal couplings of the inflaton [32] . 

Our choice of flavour charges implies that the Higgs field S decays exclusively into 
pairs of Ni neutrinos. The resulting decay rate is given by 

r°s = ^ms [l-(2M 1 /m s ) 2 ] 3/2 . (22) 

167T 
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Note that if the S boson decayed into more than one heavy neutrino flavour, this would 
lead to an interplay between lepton asymmetries which could have different signs and 
it would also change the reheating temperature. In the present paper we restrict our 
analysis to the simplest case. 



3 Entropy production through neutrino decays 

In order to understand the reheating process subsequent to false vacuum decay quantita- 
tively, we have to track the evolution of the following abundancies as functions of the scale 
factor a: the S Higgs bosons, the heavy Aq Majorana neutrinos, the B — L asymmetry, the 
standard model radiation R and the gravitinos G. The appropriate tool for this task are 
the Boltzmann equations in an expanding Friedmann-Lemaitre universe, the formalism 
of which is summarized in Appendix |Aj In what follows, we consider comoving number 
densities, 

N x (t) = a(tfn x (t), (23) 

where a is normalized to a(tpu) = 1- This quantity exhibits the advantage of being well 
defined for times prior to reheating of the universe. 



3.1 Initial conditions 

After tachyonic preheating the universe is filled by a gas of nonrelativistic S bosons as 
well as heavy A*2 and A3 neutrinos. Given the flavour structure presented in Section [2j 
the latter decay into standard model particles on time scales much shorter than the S 
boson lifetime, 

~ v 2(a- d -i) > v -2 ? i> a >0, d>\. (24) 

Hence, we do not explicitly resolve the time dependence of the ^2,3 number densities 
but approximate their evolutions by step functions, i.e. instantaneous drop-offs at times 
t 2 = t PR + 1/T% 2 and h = t PH + 1/r^ ~ t 2 (cf. Eq. g7|)), 

n Jv 2 , 3 (*) ~ q ^3 n N23 (tpH) Q(t 2 - 1) , t>t PU . (25) 

The sudden decay of the ^2,3 neutrinos sets the stage for the reheating of the universe 
which is why we choose t = t 2 — £3 as initial time when solving the Boltzmann equations. 
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Let us now determine the initial conditions at this time. When the A^ neutrinos decay 
they transfer their energy inherited from tachyonic preheating to radiation, 

Po 



Pn(t 2 



2r 



N 2 



a\t 2 



(26) 



For a thermal bath of temperature T the energy and number densities and of 
radiation quanta are given by 

^2 



7T 



PR ~ 30 g *' pT 



C(3) 



71- 



g*, n T 3 



(27a) 
(27b) 



where g* tP and g+ n denote corresponding effective sums of relativistic degrees of freedom 
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E ( T */ T ) 4 + « E ( r */ T ) 



bosons 



fcrmions 



9*,n = E ^l T f + 7 E (^/ T ) 3 ' 



(28a) 
(28b) 



bosons 



fcrmions 



In our numerical analysis we employ the values g* iP = 915/4 and g* tn = 427/2 (T; = T) for 
the minimal supersymmetric standard model (MSSM). Equating the energy densities in 
Eqs. (26) and (27a) yields the initial temperature T(t 2 ) and thereby the initial comoving 
number density Aij(£ 2 ), 

1/4 

T(t 2 ) 



30 2rjv 2 p 
n 2 g^ p a 3 (t 2 ) 

N R (t 2 ) = a 3 (t 2 )^g^ n T 3 (t 2 ). 



Note that Eq. (27) provides us with an expression for T as a function of Nr, 

7T 2 -/Vr 



T 



\C(3)0*,na s 



1/3 



(29a) 
(29b) 

(30) 



As we will argue in Section 3.2.1 this relation can be used to determine the time evolution 
of the temperature. The out-of-equilibrium decay of A^ also produces an initial B — L 
asymmetry]^] The corresponding comoving number density is given by (cf. Eq. (jl~5|)) 



N B - L (t 2 ) = e 2 N N2 (t 2 ) + e 3 N Ns (t 2 ) 
_ 2 fn u Mi 



0.2 rf 



-N N2 (t 



(31) 



'EW 



l ln Ref. |28j only the initial B — L asymmetry from tachyonic preheating is taken into account. 
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where the comoving number density N^ 2 at tpH follows from Eq. (20b). 

At t = t 2 the dominant contribution to the energy density resides in the gas of non- 



relativistic S bosons (cf. Eq. (20a)), 



Ps{h 



1 



a 3 (t 2 



(1 - 2r N2 )p . 



which corresponds to an initial comoving number density 

Ps(t 2 ) 



N s (t 2 ) = a s (t 2 )n s (t 2 ) = a s (t 2 ) 



m s 



(32) 



(33) 



According to Eq. flA.7 ), the initial phase space distribution function fs(h,p) can be 
inferred from ns(t 2 ). Guided by the results of Ref. [28] we make the ansatz of a delta- 
peaked momentum distribution, i.e. fs(t 2 ,p) oc S(p), which leads to 

5(k) 



f s (t 2j p) = 2ir 2 N s (h 



k' 2 



k = a(t 2 )p(t 2 ) = a(t)p(t) . 



(34) 



Because of the chosen mass hierarchy Mi <C ms — M 2> 3, the amount of N\ neutrinos 
produced during tachyonic preheating as well as through S decays up to t 2 is negligibly 
small. Likewise, in absence of standard model radiation no gravitinos are produced until 

h, 



N Nl (h 







0. 



(35) 



The time dependence of the scale factor a(t) is governed by the Friedmann equation. 
For a flat universe and constant equation of state u = p/p between some time t and t, 
one has 



a(t) = a (i( 



UJ 



8tt 



1/2 



Ptot (*o) (t-t ) 



3(l+w) 



(36) 



After preheating, until the decay of the S bosons around t$ = t 2 + l/T° s , the system is 
mostly matter-dominated. In view of the mass hierarchy Mi/mg ~ rf 1 the subsequent 
S decay into relativistic N\ neutrinos, however, entails a continuously changing equation 
of state. We account for that behaviour by working with two constant effective equation 
of state coefficients ui 2 and us for times t 2 < t < t$ and t > ts, respectively. For times 
£ph < t < t 2 , we take u — 0. u 2 can be deduced from the decrease in the total energy 
density until t$. Keeping only the leading contributions to p tot one has 

3(l+w 2 ) 

(37) 



po/a(t 2 ) 3 



a(t 2 



12 



With explicit expressions for p$ and pn 1 at hand (cf. Sections 3.2.2 and 3.2.3) this 
equation can be solved numerically for w 2 - Within the region in parameter space to which 
we restrict our study (cf. Section [2]), 002 typically turns out to be closer to than to 1/3 
indicating matter domination. The u 2 mean value, standard deviation, minimum and 
maximum values are found to be 

u 2 = (2.4 ± 1.6) x HT 2 , uf n ~ 8.8 x 1(T 3 , u™ 3 * ~ 8.7 x 10~ 2 . (38) 

At times t > t,s most of the initial energy density has already been transferred to Ni 
neutrinos or subsequently to standard model radiation and we can safely use u)g — 1/3. 

3.2 Boltzmann equations 

Using the conventions introduced in Appendix [A] we now write down the Boltzmann 



equations relevant for our specific scenario. In Section 3.2.1 the Boltzmann equations for 



all species but the gravitino (S, Ni, B — L, R) are given for phase space distribution 



functions. Subsequently, they are discussed one by one in Sections |3.2.2| through |3.2.5 
For the gravitino component we directly give the integrated Boltzmann equation for the 
comoving number density in Section 3.2.6| 



3.2.1 Collision operators 

The dynamics of our system is dominated by three types of particle interactions: S boson 
decays into pairs of N\ neutrinos, N\ neutrino interactions with standard model lepton- 
Higgs pairs £H and £H and supersymmetric QCD 2—7-2 scatterings responsible for the 
production of gravitinosj^] 

The Boltzmann equations for the S bosons and the Ni neutrinos take the form 

Cfs = C s {3^N 1 N 1 ), (39) 
Cf Nl = 2 C Nl (S NxNi) + C Nl (Ni £H, £H) , (40) 



where the factor of 2 in Eq. (40 ) accounts for the fact that two N% neutrinos are created per 



S decay, and Cn 1 (Ni -h- £H,£H) encompasses Ni decays into particles and antiparticles, 

C Nl (N x ^ £H, £H) = C Nl (Ni ^ £H) + C Nl (JVi ^ IS) . (41) 

Since we expect their effects to yield only minor corrections, we do not include the rescat- 
terings of N\ neutrinos into S bosons as well as scatterings involving massive Z' bosons 



5 For notation convenience, we refer to H u as H from now on. 
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which arise in the course of B — L breaking. We have checked that the scatterings 
Nil -H- qu, Niu -H- £q and N%q -H- £u do not affect the final B — L asymmetry. 

The Boltzmann equation for the B — L asymmetry is defined in terms of the respective 
equations for lepton number L and anti-lepton number L 

Cf B -L = C(f L -f z ) , (42a) 
Cf L = C L (£H Aq ) + 2 C r L cd (IH £#) , (42b) 
£/ £ = C L (£H Nx) + 2 C| cd (Zff o . (42c) 

The collision operators for decays and inverse decays are able to mimic AL = 2 scatterings 
of the type £H ^ Aq ^ £H with on-shell N\ neutrinos in the s-channel. They, however, 
ignore off-shell scatterings even though these will equally affect the final asymmetry. This 



leads us to adding reduced collision operators C£ cd and C T £ d to Eqs. (|42bJ) and (|42cj) that 



account for the production and decay of off-shell neutrinos, £H ^ N* ^ £H. In Ref. 
it has been shown that the on- and off-shell contributions to the total CP asymmetry in 
£H -H- £H scatterings cancel up to O ((h^) ). Hence, one may equivalently say that the 
reduced collision operators subtract scatterings with real intermediate states. 

The temperature T of the thermal bath can be determined as a function of the scale 
factor from the covariant energy conservation 

p tot + 3H (p tot + p tot ) = . (43) 

After inserting the explicit expressions for the total energy and number densities this 
relation becomes a non-linear first-order differential equation for T. For simplicity, we 
assume that the energy transfer to the thermal bath happens instantaneously, p tot ~ 0. 
Consider a spatial volume V in which Aq neutrinos of average energy decay into 
lepton-Higgs pairs. Per decay the energy density of the thermal bath is then increased 
by s^x/V and a new thermal equilibrium at a slightly higher temperature is established 



right after the decay. According to Eq. (27) the latter entails an increase in tir, 



( \ , £ Ni i ^ ^TVi . r R /AA x 

nR ^ nR \ l + W*J ^ nR+ 4V^r R = nR+ v- (44) 

Therefore, producing two standard model particles adds r R = 3^/ (4p fl /n R ) radiation 
quanta per unit volume V to the thermal bath. This leads us to an effective Boltzmann 
equation for the number density of radiation 

Cf R = r R (C e (£H N x ) + C- t {£H N x )) , (45) 
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where all CP violating contributions have been neglected. 

Having formulated the Boltzmann equations for distribution functions we now calcu- 
late the various collision operators explicitly and simplify the equations as far as possible. 



3.2.2 Higgs bosons S 

In terms of the S decay rate 

T ° s = 2^~sl dU{Nl > Nl) ( 27r )^ (4) (&out - 5>m) \M (S iV^OI 2 , (46) 
the collision operator for S decay is given by 

(47) 

This is a linear homogeneous ordinary differential equation which has a unique solution 



Cfs = C S (S -> N.N,) = ~T° s f s . 



for each initial value. Given the initial distribution function fs(t2,p) m Eq. (34) we find 



fs{t,p) = / s (i 2 ,p)exp 



(48) 



-m s Y% j dt'E s \t' 

*2 

E s {t') = ^{a/a'fp 2 + m| , a = a{t) , a' = a(t') . (49) 
Thanks to the momentum delta function in fs(ti,p) the time integration becomes trivial, 
fs(t, P ) = fs(h,p) e-^W = 2^N s {t 2 ) 5 ^e-r°s^) , k = a (t)p. (50) 



Hence, the comoving number density N$ simply falls off exponentially, 

N s (t) = N s (t 2 )e- r ^-^ . 

3.2.3 Heavy Majorana neutrinos Ni 



(51) 



The Boltzmann equation for Ni neutrinos also involves a collision operator for S decay. 
Using the tree-level amplitude squared 



\M(S -> NxNi)\ 2 = 2(/i?) 2 m| [l - (2M 1 /m s ) 2 ] 



(52) 



the general collision operator in Eq. (A. 3) takes the form 

,n\2^2 



C Nl {S^N 1 N l ) = ^-^[l-{2M 1 /m s f] [ dtl^l^; S)frr)W> f s . (53) 

J 
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With the explicit expression for fg in Eq. (50) the phase space integration becomes 

1 ttN s 



cfn(Aq|Aq;S) (2vr) 4 ^/ s 



-5(E Nl -m s /2) 



(54) 



a 3 8m s E Nl 

where a symmetry factor 1/2 results from the fact that two Aq neutrinos are involved in 
the decay process. Employing the result Eq. (22) for the S decay width we obtain for the 
collision operator 



C Nl (S -> A^Aq 



- 3 ^f^ [l ~ (2M 1 /m s ) 2 ]- 1/2 HE Nl -m s /2) . (55) 



The collision operator for Aq decay into standard model particles has the familiar form 



C Nl (Aq tH, £H) 



Mi 
'E~ 



Ni 



(56) 



In total, the Boltzmann equation for Aq neutrinos encompasses two production and one 
decay term. On the one hand, the collision operator in Eq. (55) and the term proportional 
to fff in Eq. (56) represent Ni production from S decay and from inverse decays in 
the thermal bath, respectively. On the other hand, the term proportional to in 
Eq. (56) accounts for Aq decays. It is convenient to decompose the Aq population into 
two independently evolving components: Nonthermal neutrinos Aq stemming from S 
decays and thermal neutrinos N± originating from the thermal bath. The respective 
Boltzmann equations are then given as 



Ni 



T Mi eq 

1 Ni J Ni ' p 1 NiJNi ■ 



E N 



' X ' " [l-(2M 1 /m s ) 2 Y 1/2 5(E Nl -m s /2), (57a) 



Mi 
Em 



(57b) 



The sum of the nonthermal and thermal distribution functions yields the total Aq distri- 
bution function, f Nl = f Ni + 

The Boltzmann equation for the nonthermal neutrinos Aq 5 can be solved exactly. 
Starting from zero initial abundance we find 



f Nl (t, p) = 2n 2 T s [1 - (2M 1 /m s ) 2 ] ^ J dt' 



t-2 



8(E Nl (t')-m s /2) 



x 



J^^[-M,T%JdfE N l { t 



(5f 
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where E Nl (t') is defined analogously to E s (t') in Eq. (49). The energies E Nl therefore 
redshift as 



Em (t) = E Nl (t 



a 



1 + 



Mi 



1/2 



(59) 



If we evaluate this relation with E^^t') = ms/2 we obtain the redshifted energy at time t 
of a neutrino N± that has been produced in S decay at time t'. Let us denote this quantity 
by S Nl (t',t), 



£nx (t , t) 



ms a 



1 + 



2M X 

m s 



1 1/2 



(60) 



The energy delta function in the integrand of Eq. (58) thus turns E Nl (t") into S^^t', t"). 
Meanwhile, it can be rewritten as a function of E Nl (t), the energy at time t, 



5{E Nl (t')-m s /2) 



a'\ 2 777,5/2 



a) £ Nl (f,t) 



5(E Nl (t)-S Nl (t',t)) 



(61) 



The final result for then reads 

t 

f Ni (t,p) = l 2 vr 2 r° [{m s /2f - Ml]' 112 j dt> 

x ^f^y ex P {-M x v\j dt"£- N l{t',t") 



-S(E Nl (t)-£ Nl (t',t)) 



(62) 



By integrating over the Nf phase space, we obtain the number density nff of the non- 
thermal neutrinos 



d 3 P fS _ 2r° 5 



x ) (27T) 



a* 



(63a) 



J dt' Ns^expl-M^J dt"S N l(t\t") 

dt'5n s Nl (t',t) . (63b) 
The corresponding results for the energy and interaction densities pfj and 7^ = ^(Nf — > 



t-2 



17 



£H,£H) can conveniently be expressed using Snff (t, t') introduced in Eq. (63b), 

t 

p s Nl = [ dt'£ Nl (t\t)5n s Ni (t',t) } (64) 



*2 

t 

ik = J dt ' SN M ( t',t) T ° N ^ t] = n ^ ' (65) 

where denotes the N\ decay width weighted with the average inverse time dilatation 
factor for nonthermal neutrinos 

r - =<H>/». =i/ rft, OT^ (t ' t)r °- (66) 

The exact phase space distribution function /jy for thermal neutrinos Nf is given as 
the unique solution of Eq. (57b) for the initial distribution (t2,p) = 0, 

fl(t, P ) = J dt'exp (—MiT° Nl j dt"E- N l{t") J T%f%tf,p) . (67) 

As the thermal neutrinos are produced within a broad range of energies, it cannot be inte- 
grated over phase space as simply as in the nonthermal case. However, since the thermal 
neutrinos inherit their momentum distribution from the thermal bath it is reasonable to 
assume that they are approximately in kinetic equilibrium, 

fl(t,p) « T$/£(f,p), m,p)=e~ E ^ T . (68) 



This approximation holds if the quotient f^/f^, with taken from Eq. (67), is inde- 
pendent of the neutrino momentum p. In Appendix [B] we will demonstrate numerically 
in the context of a specific parameter example that this can usually be assumed to be the 
case in our scenario. Under the assumption of kinetic equilibrium the comoving number 
density Nj^ is the unique solution of the integrated Boltzmann equation 

aH ^ N l = - K - Affl r& • (69) 

Here, stands for the N% decay width weighted with the average inverse time dilatation 
factor for thermal neutrinos 
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where z = Mi/T, and K\p{z) are modified Bessel functions of the second kind. Note that 



7 ^ = 7 (JVf -> £ff, Iff) = n^I% , (71a) 

7 £ = 7 eq (#i Iff) = n^r* . (71b) 

3.2.4 B — L asymmetry 

The collision operators for decays and inverse decays present themselves as 

C L (Iff ~ l-M (JVx £ff)| 2 - l-M (Iff iVOI 2 , (72a) 

C z (Iff iVx) ~ f Ni \M (Ni -> Iff) | 2 - |M (Iff N t ) | 2 . (72b) 

Using the definition of the CP parameter ei and CPT invariance, the various partial 
amplitudes squared are related to the total amplitude squared as follows 

\M (N x -> £H)\ 2 = \M (Iff ^N l )\ 2 = l -{1 + ex) \M Nl | 2 , (73a) 

\M {N, -» Iff) | 2 = \M {£H -> iV\)| 2 = X - (1 - ei) lA^J 3 , (73b) 

where, at tree-level, |A^jvJ 2 is given as 

\M Nl \ 2 = \M (Ni -> £H,£H)\ 2 = 4(^ t ^) n Mj 2 . (74) 



The reduced collision operators in Eq. (42) account for the production of off-shell 



neutrinos N± which subsequently decay into the CP conjugate of the lepton-Higgs pair 
from which they were produced. Working up to leading order in e\ we may take the 
decays AT* — > £H, £H to equally branch into particles and antiparticles 



C T L ed L {IE <-> £H ) ~ ± 



fifa ■ \ \M (Iff -> N*)\ 2 - Mh ■ I \M (Iff -> iV- 



(75) 



For M x <C 10 14 GeV the CP preserving parts of the off-shell scatterings are negligibly small 
[34] . We thus discard them keeping only the CP violating contributions to the reduced 
collision operators. Imposing that the total CP asymmetry of lepton-Higgs scatterings be 
zero up to O ((/^i) 4 ), we deduce 

\M{£H ^N{)\ 2 = - l -{-t l )\M Nl \ 2 , \M{£H Nl)\ 2 = - l -{+e l )\M Nl \ 2 . (76) 
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The above results allow us to write all collision operators as integrals over the total 
amplitude squared I-MatJ 2 . Assuming kinetic equilibrium for leptons and antileptons as 
well as thermal equilibrium for all other standard model particles, we obtain 



tf B -L = 7^~ [ dUiilH;^) (2vr) 4 ^ \M Nl \ 



(f Nl - a) - 



Nb-l f ,,i 
Nf q ' 



' f eq 



■ (77) 



As for the Aq Boltzmann equation, we split the Aq distribution function into its thermal 
and nonthermal parts. After integrating over phase space we retrieve the interaction 



densities 7^, 7^ and 7^ (cf. Eq. QA.8[ )) 



Q-H-^-Nb-l = a 3 
da 



N, 



2n: 



B-L eq 
eq 



(7E 



In terms of comoving number densities and averaged decay rates this Boltzmann equation 
then reads 



vH^Nb-l = e l N^ l T s Nl + e 1 {N T Ni - N%) T T Ni - ^ ^Nb-l ■ 



2A, cq 



(79) 



Similarly to the Aq abundance we also split the B — L asymmetry into two components: 
A nonthermal asymmetry N§_ L produced in Af decays and a thermal asymmetry Aj_ L 
generated from the thermal bath, 



d 



aH da N B-L-^N^T s Nl 2.\T" ' v '^ ; - i ' 



N eq 

JV ^i pT N S 



d 



aH Ta N T B _ L = e, (A^-A^Jr 



N' 



eq 

Ni T T N T 



T 

Ni 2N eq A x 1 1 L1 ~ L 



(80a) 
(80b) 



A comparison of the corresponding final baryon asymmetries 77^ and 77^ will allow us 
to identify the relative importance of nonthermal and thermal leptogenesis in different 
regions of the parameter space (cf. Section [5]). 



3.2.5 Radiation R 

In order to obtain an effective Boltzmann equation for the number density of radiation 
quanta we add up the contributions coming from decays into standard model particles 



and antiparticles. Neglecting any CP violating effects Eq. (45) can be written as 



tfR = TR ^p J rfn W^ Ar i)(2^) 4 ^ (4) ^l^ 1 | 2 [2^ 1 -2/a ■ (81) 
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Splitting into its thermal and nonthermal parts and integrating over phase space, one 
gets 



d 



aH Ta NR = r * JV ^ I % + r « - ^) r ^ ' 



52) 



Since the two different sorts of Ni neutrinos possess different average energies, two inde- 
pendent factors and have been introduced in order to keep track of the radiation 
quanta produced in the decays of nonthermal and thermal neutrinos respectively, 



3f T 
T _ 0t N! 



R ~ ±SR ' 



R 



1 



R 



The average energies per (non-)thermal neutrino as well as the average energy per radia- 
tion quantum are obtained from the respective ratios of energy and number densities (cf. 



Eqs. (27), (63) and (64)), 



£r = PR/n R , = p S Nl /n S Nl , e T Nl = p T N Jn T Nl = 3T + 



K 2 (z) 



Mi. 



i4) 



rj| and r\ clearly depend on the temperature T which, in turn, is deduced from the 



radiation number density N R according to Eq. (30). 



3.2.6 Gravitinos G 

Gravitinos are produced through scattering processes in the thermal bath. The evolution 
of their comoving number density is governed by the Boltzmann equation 



S5) 



The dominant contribution to 7g comes from QCD scatterings. For supersymmetric QCD 
and up to leading order in the strong gauge coupling g s , one has [15] 



1g(T) 



1 + 



m 



J(T)\ 54C(3)^ 2 (T) 



3m~ 



tt 2 M 2 



rpb 



In 



T' 2 



m 2 (T) 



+ 0.8846 



36) 



Here, nig denotes the energy scale-dependent gluino mass and m„ is the gluon plasma 



mass, 



m- g {T) 



9 2 S (T) 
9s (Po) 



m s (jio) 



in, 



,(T) = \fm g s (T)T. 



(87) 
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Point 



Label v B ^ L [GeV] M x [GeV] m x [eV] 



m 5 [GeV] 



m? 



[GeV] 



Section 


i] Red circle 


5.8 x 10 13 


1.4 x 10 10 


3 x 10~ 3 


100 


800 


Ref. 


White circle 


3.0 x 10 12 


1.0 x 10 10 


1 x 10~ 3 


100 


800 



Table 2: Values of the input parameters chosen for the discussion in Section [i] and Ref. \S$, respectively. 
In Figs. [5[ [6| @ [T3| and 14 in Section [5] and Appendices and [D] tie positions of both parameter 
points in parameter space are marked with different labels as indicated. 



where we choose the Z boson mass M z as reference scale /i . The scale dependence of g s 
is dictated by the corresponding MSSM renormalization group equation 

-1/2 



9s{p) = 9s (//c 



3 2 
1 + 5T 9s (//o) In (Wo) 



57T 2 



with /x being the typical energy scale during reheating. It can be estimated by the average 
energy per relativistic particle in the bath: /i ~ Er ~ 3T. For instance, at temperatures 
T = 10 8 , 10 10 , 10 12 GeV the strong gauge coupling takes on values g s = 0.90,0.84,0.80. 
The gravitino mass mg and the low-scale gluino mass rrig (/i ) remain as free parameters. 

4 An illustrative example 

In the previous section we have seen how the decay of the false vacuum of unbroken 
B — L symmetry generates - via production and decay of heavy neutrinos - entropy, 
baryon asymmetry and dark matter. We have numerically solved this inital- value problem 



by means of Boltzmann equations, with the initial conditions described in Section 3.1 
Before we turn to a detailed discussion of the parameter space we first describe, as an 
example, one solution for a representative choice of parameters. A similar study, albeit not 
as detailed, has already been performed in Ref. [3J. The values for the input parameters 
chosen in this section as well as in Ref. [3J are outlined in Tab. [2] In the present case, 
the selected values for vb-l and M\ correspond to Froggatt-Nielsen charges of a = 1/2 
and d ~ 1.5. Note that we have adjusted M\ such that we obtain the right gravitino 
abundance for dark matter. The input parameters in Tab. [2] directly determine a couple 
of further important parameters 

m s , M 2 , 3 ~ 4.1 x 10 12 GeV , A ~ 5.0 x 10~ 3 , e x ~ 1.4 x 10~ 6 , e 2 , 3 ^ -4.1 x 10~ 4 . (89) 



We have chosen opposite signs for the CP asymmetries e\ and €2,3, so that one can easily 
distinguish their respective contributions to the final B — L asymmetry. 

22 



Inverse temperature M\/T 
10" 1 10° 10 1 10 2 




Scale factor a 



Inverse temperature M\/T 
10" 1 10° 10 1 10 2 




Scale factor a 



Figure 1: Comoving number (upper panel) and energy (lower panel) densities for S bosons, Ni 
neutrinos produced in S decays (Nf), thermally produced Ni neutrinos (N±), N\ neutrinos in thermal 
equilibrium (N° q , for comparison), radiation (R), B — L charge, and gravitinos (G) as functions of the 
scale factor a. The corresponding values of the model parameters are given in Tab. [2[ 
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4.1 Comoving number and energy densities 

The evolution with the scale factor a of the comoving number densities and the various 
components of the energy density are presented in Fig. [T} Both plots start at the time of 
^2.3 decay, a (£2) — 1-04, and end at a final scale factor of a/ ~ 2.53 x 10 5 . The values of 
the scale factor corresponding to S and N\ decays as well as to reheating (cf. Appendix |C| 
are also indicated 

a (t s = t 2 + l/r°) ~ 400 , a (t, = t s + l/T° Nl ) ~ 1200 , a (t RH ) ~ 3000 . (90) 

Tachyonic preheating results in an initial state ai t = t 2 that mainly consists of S 
bosons and to a smaller degree standard model radiation that stems from the decay of 
the ^2,3 neutrinos, pn (£2) / Ps (^2) ~ 10~ 5 , and which inherits an initial B — L asymmetry 
equivalent to T]b — —3.9 x 10~ 5 . Around t = ts the S bosons decay into relativistic 
and nonthermal N± neutrinos. Their subsequent decay into standard model particles 
then washes out the initial (negative) asymmetry, builds up a new (positive) asymmetry 
and leads to the production of the main part of the radiation. The energy transfer to the 
thermal bath, i. e. the process of reheating, takes place between a l RH ~ 27 and a^ RH ~ 5800. 
At these values of the scale factor the derivative of the comoving radiation energy density 
a 3 pn vanishes. Meanwhile, thermal neutrinos and gravitinos are continuously produced in 
the thermal bath. As both species inherit their momentum distributions from the bath, 
they are always in approximate kinetic equilibrium (cf. Appendix [B]) . From a ~ 4100 
onwards the number density of the thermal iVi neutrinos exceeds the thermal equilibrium 
abundance. At a ~ 10 5 the B — L asymmetry and the gravitino abundance have reached 
at their final values. 



4.2 Baryon asymmetry 

The present value of the baryon asymmetry as well as its nonthermal and thermal con- 
tributions are obtained from 



r>° O N S,T 

n B S , T S,T y*,s iv B-L 

Vb = -q =Vb + Vbi Vb = c s P h 



(91) 



a f 

In the supersymmetric standard model the sphaleron conversion factor is c sp h = 8/23, the 
effective number of degrees of freedom at high and low temperatures is g* jS = 915/4 and 
g® s = 43/11, respectively, and the number density of photons is iV 7 = a 3 g 7 £(3)/7r 2 T 3 . 
For our choice of parameters we obtain the asymmetries 

r] B ~ 1.9 x HT 8 , r?f ~ 1.9 x 10" 8 , r} T B ~ 2.8 x 10~ 10 . (92) 
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The calculated baryon asymmetry is larger than the observed one, r] B bs ~ 6.2 x 1CT 10 [20] . 
by about a factor 30. This is consistent since e\ is an estimate for the maximal CP 
asymmetry. We find that r/s is dominated by the nonthermal contribution due to iV-f 
decays, vfe ~ rjB- The contribution from thermal neutrinos, even though it reaches the 
right order of magnitude, slightly falls short of the observed value. 

Let us emphasize that given the choice of parameters in Tab. [2] standard thermal 
leptogenesis, with a given thermal bath, is able to produce the right amount of baryon 
asymmetry. Using a final efficiency factor of Kf(rh\) ~ 0.1, one obtains (cf. [33J) 



3 0° 

Vb = 7— Csph ei «/(mi) ~ 6 x 10- 10 . (93) 

In the case under study rj B turns out to be smaller than rf B roughly by a factor of 2 because 
the entropy production during iVf decay enhances the washout rate due to inverse N\ 
decays. 

The evolution of the nonthermal and thermal lepton asymmetries is controlled by three 
different interaction rates which enter the Boltzmann equations (80a) and (80b) 

N s N T - N eq N eq 

1 B-L - e l J^s 1 ATi > 1 B-L - e l jyr 1 ^1 ' 1 W ~ 1 N x ■ l y4 J 

_B — L B — L £ 

They account for the decay of nonthermal neutrinos, the decay of thermal neutrinos and 
the washout effects due to inverse neutrino decays, respectively. Their relative importance 
as well as their influence on the generation of the asymmetries are illustrated in Fig. [2j 
The respective interactions become efficient once the corresponding rates are of the same 
order as the Hubble rate or larger. This is why it takes until ai ~ 210, when Tg_ L /H > 
(9(1) for the first time, for A?J_ L to begin to increase. At 02 — 480 the initial negative 
asymmetry has been compensated by the generated positive one, and A?J_ L changes sign. 
Subsequently, for scale factors around a(ti), the ratio V B _ L /H remains approximately 
constant leading to the generation of the main part of the asymmetry. Meanwhile, due 
to the continuous entropy production from nonthermal neutrino decays, the washout 
processes gain in importance. At a 3 ~ 1500 the rate Ty/ becomes comparable to H, which 
is reflected in a slight decrease of the slope of N^_ L . From a 7 ~ 6300 onwards, which is 
shortly after p# = pf fi , the washout even dominates over the asymmetry production from 
Ni decays. Hence, the maximal nonthermal asymmetry reached at aj is slightly washed 
out until it eventually freezes out when Tw drops below H at axo — 10000. Notice that 
^b-l aheady becomes irrelevant at ag ~ 8300. 

The decays and inverse decays of thermal neutrinos lead to a continuous production 
of a thermal asymmetry with a negative sign until the rate Tg_ L is overcome by Ty/ at 
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Fi gure 2: Evolution with the scale factor a of (upper panel) the interaction rates r^_^ J and 

Tw normalized to the Hubble rate H, and of (lower panel) the nonthermal and thermal parts N§_ L 



and Ng_ L of the generated B — L asymmetry. The rates are defined in Eq. (94), and the asymmetries 



were introduced in Eqs. ( 80a ) and ( 80b I . The dashed lines and the integer numbers above the top frame 
edge in the lower panel refer to the various values of the scale factor and their numbering as used in the 
discussion of this figure in Section \4.2\ The gray band in the upper panel indicates where the interaction 
rates are of the same order as the Hubble rate H . 
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a4 — 3500. Following that moment, washout processes push N R _ L back to N R _ L = 0. 
This development is reinforced by thermal neutrino decays once iVjy has exceeded the 
equilibrium number density at as ~ 4100. Until a§ ~ 6500 the thermal asymmetry 
is then converted into a positive asymmetry. After a§ ~ 4700 the rate T B _ L permanently 
dominates over rv, and the thermal asymmetry does not decrease anymore after as- 
Instead it freezes out at its maximum value when Y T B _ L jH ~ (9(0.1) at an ~ 12000. 

Finally, parallel to the onset of entropy production at a RH ~ 27, the rate T B _ L exceeds 
Tw, and T B _ L slightly increases. 



4.3 Radiation temperature 



Having solved the Boltzmann equation (82) for the number density of radiation quanta, 



we obtain from Eq. (30) the evolution of the plasma temperature T which is plotted 
in Fig. [3] We find that the reheating process between a RH ~ 27 and a RH ~ 5800 is 
accompanied by an approximate temperature plateau around T ~ 6 x 10 9 GeV. Especially 
until S boson decay around a{t$) — 400 the temperature is essentially constant. This 
is due to the continuous production of nonthermal neutrinos which do not efficiently 
decay before a (ti) ~ 1200. With nonrelativistic S bosons still representing the dominant 



contribution to the energy density, the comoving N^- number density approximately scales 



like oc f t dt' oc a 3 / 2 . According to the Boltzmann equation (82) for radiation, the 



comoving number density Nr then grows like the volume, implying a constant temperature 
aH^-N R oc oc a 3/2 , N R oc a 3 , T = const. (95) 

Once the production of nonthermal neutrinos ceases, not as much radiation is produced 

anymore and the temperature begins to drop. During the phases of adiabatic expansion 

T decreases like the inverse of the scale factor, T oc 1/a. 

The actual reheating temperature Trh is reached once the Hubble rate H becomes as 

s 

Ni 



small as the effective decay rate Tfj of the nonthermal neutrinos (cf. Appendix 



1% (tjiff) = H(t RH ) , T RH = T(t RH ) . (96) 

For the chosen set of parameters this happens at a (tRjj) ~ 3000, with H = Tfj = 52 GeV, 
and the corresponding temperature turns out to be 



RH 



4.1xl0 9 GeV. (97) 



A detailed discussion of how this result for the reheating temperature can be estimated 
on the basis of the input parameters is given in Appendix [Cj 
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Figure 3: (Upper panel) Temperature T and (Lower panel) the nonthermal and thermal correction 
factors and as functions of the scale factor. The temperature is calculated from Eq. (30 1, the 



correction factors were introduced in Eq. (83) and are used in the Boltzmann equation (82) for radiation 



to ensure a correct counting of radiation quanta. 
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The lower panel in Fig. [^presents the evolution of r|| and as functions of the scale 
factor a, the two correction factors which effectively keep track of the average energy 
per nonthermal / thermal neutrino in relation to the typical radiation energy er 



as discussed in Sections 3.2.1 and 3.2.5 Until the onset of reheating, steeply rises. 
This is the consequence of an adiabatically dropping temperature and, on top of that, 
the progressively increasing effectiveness of the S boson decays which push the average 
Nff energy more and more towards mg/2. Between a % RH and a(ts) the temperature stays 
rather constant and, as we have checked numerically, eff has saturated close to ms/3. 
Hence, only varies little around a value of ~ 50 during that time. After a(ts) the S 
boson decays become less frequent, the energies are redshifted and decreases. This 
trend is stopped around a(£i) when the temperature begins to fall again and the decay 
of the nonthermal neutrinos themselves sets in. These decays tend to remove rather 
long-lived and hence redshifted neutrinos from the spectrum leading to an increase in 
eff . Finally, after reheating the evolution of is again dominated by the adiabatically 
decreasing temperature. 

The initial value of the thermal correction factor, r^(£ 2 ) — 0.78, is close to 49/60 and 
hence to what is expected for a relativistic fermion coupled to the massless degrees of free- 
dom of the MSSM. The fact that it is even a bit larger is due to the negligible imprecision 



of calculating by means of classical statistics (cf. Eq. (84)). Once the temperature has 



dropped below M\ around a ~ 17, the thermal neutrinos become nonrelativistic and 
increases. This evolution is only delayed by the constant temperature during reheating. 
After reheating continuous to increase like oc a. 

4.4 Gravitino dark matter 

The present contribution from gravitinos to the total energy density is given by 



Pc g*, s A7 



G 



(9* 



where p c = 1.052 x 10~ 5 h 2 GeVcm -3 is the critical density and n° = 410 cm -3 is the 
number density of CMB photons. Recall that after fixing all other parameters we have 
chosen M\ such the gravitino abundance equals the observed one for dark matter 

fig/1 2 ~ 0.11. (99) 

In Appendix [D] we demonstrate that this result can be easily reproduced by means of 
certain semi-analytical estimations. 
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Figure 4: Interaction rates = Jq/tIq, 
radiation, normalized to the Hubble rate H. The rates r|j and were introduced in Eq. (100). The 
gray band indicates where the interaction rates are of the same order as the Hubble rate H. 



The evolution of the gravitino abundance with time is controlled by the production 



rate = 7g/^g (cf. Eqs. (85) and (p6J)), which strongly depends on the temperature. To 
demonstrate the close relation between the production of gravitinos and that of radiation 
we plot in Fig. together with the rates and at which radiation is produced 



due to nonthermal and thermal neutrino decay 



R R at Ni ■> L R~ r R , r N\ ' l iuu J 



As expected we find that only the decay of the nonthermal neutrinos efficiently influences 
the radiation abundance. Between a ~ 22 and a ~ 7400, which is basically identical to 
the time interval in which reheating takes place, is of the same order as the Hubble 
rate. T S R /H reaches its maximal value before S boson decay around a(tg), subsequently 
decreases a bit and finally drops off shortly after a RH . The ratio T^/H is at most of 
order O(10~ 2 ) which is the case towards the end of reheating when the thermal neutrinos 
are close to thermal equilibrium. The rate traces the efficiency of the nonthermal 
neutrino decays: While the iVf decays are not active yet, Fq/H decreases due to the 
falling temperature. But as soon as rfj becomes competitive with the Hubble rate around 
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a l RH , Tq/H bends over and eventually it reaches values of order 0(1) and even larger. On 
the other hand, once the nonthermal neutrino decay has ended, returns to its ordinary 
behaviour that we expect for adiabatic expansion. In total gravitino production occurs 
between a ~ 79 and a ~ 18000. The main part of the gravitino abundance is, hence, 
produced towards the end or after reheating. 



5 Results and discussion 



The parameter point selected in the previous section was chosen such that we readily 
obtained the right baryon asymmetry and gravitino abundance. Now we extend our dis- 
cussion to a quantitative analysis of the entire parameter space and determine the bounds 
within which consistency between successful leptogenesis and gravitino dark matter can 
be reached. According to the flavour model introduced in Section [2] we are free to vary 
the neutrino mass parameters Mi, rh\ and vb-l- On the supergravity side the gravitino 



and gluino masses mg and rrig represent free parameters (cf. Section 3.2.6). Moving in 
parameter space changes the interaction rates relevant to our scenario, most notably the 
production and decay rates of the N\ neutrino. This has consequences for the reheating 



process (cf. Section 5.1 ), the generation of the baryon asymmetry (cf. Section 5.2 ) and the 



thermal production of gravitinos (cf. Section 5.3). By imposing the two conditions 



VB = VB+Vl > V°B S ~ 



6.2 x 10~ 10 



Q 5 h 2 



Q DM h 2 ~ 0.11 



(101a) 
(101b) 



we are able to identify the regions in parameter space in which both, the present baryon- 
to-photon ratio and the dark matter density are successfully generated. In this manner, 
we obtain a link between neutrino and superparticle masses. The parameter dependence 
of the reheating temperature and the interplay of nonthermal and thermal leptogenesis 
follow along the way. 



From the allowed range for the B — L breaking scale (cf. Eq. (10)), we consider the 



boundary values and an intermediate scale. All three values are associated with different 
ranges for the heavy Majorana mass Mi (cf. Eq. (J7a))), 



VB _ L = 3.4 x 10 12 GeV 
v B -L = 5.8 x 10 13 GeV 
v b _l = 1.0 x 10 15 GeV 



1.3 x 10 5 GeV < Mi < 1.1 x 10 10 GeV, 
2.1 x 10 6 GeV < Mi < 1.9 x 10 n GeV, 
3.7 x 10 7 GeV < M 1 < 3.3 x 10 12 GeV . 



(102a) 
(102b) 
(102c) 
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Furthermore, in order to take into account the 0(1) uncertainties in the Yukawa couplings 
h v , we allow the effective neutrino mass rh\ to vary in the range 

l(T 5 eV < mi < 10 -1 eV. (103) 

The gravitino mass is taken from the interval 

30 MeV <m d < 700 GeV . (104) 

In view of the present bound on the gluino mass, m g > 700 GeV, imposed by collider 
searches [351136] , we use a mass of m g = 800 GeV as a representative value in this section. 
Different choices of m g would lead to similar qualitative results, the only difference being 



that all values of mg would have to be rescaled (cf. Section 5.3 and Appendix D). 

In all plots of the parameter space presented in this section as well as in the Appen- 
dices [C] and D (Figs. [5} [6} [7} [8} [ToJ [I3] and 14) we mark the positions of the two parameter 



points listed in Tab. [2] The point the discussion in Section [| was based on as well as the 
point that was investigated in Ref. [3]. 

5.1 Reheating temperature 

The concept of temperature is only applicable as long as the interactions in the system 
under study are in thermal or, at least, kinetic equilibrium. Hence, regarding our scenario, 
it is not before the creation of an initial thermal bath due to quick thermalization of the 
^2,3 decay products that we can meaningfully speak about a temperature. Subsequently, 
the main part of the energy density continues to reside in nonthermal particles. At first 
most of the energy density is carried by the S bosons, and then, from t — is onwards, by 
the nonthermal Aq neutrinos. The energy transfer to the thermal bath, i. e. the reheating 
of the universe, becomes fully efficient when the Aq neutrinos decay into standard model 
radiation. In first approximation, this happens once the Hubble rate H has dropped to 
the value of the effective decay rate of the nonthermal Aq neutrinos, 

1% (tun) = H(t RH ) , T RH = T(t RH ) , (105) 

where r^- is the vacuum decay rate weighted with the inverse time dilatation factor 



as defined in Eq. (66). The reheating temperature Trh can then be obtained by applying 
Eq. ( 105 ) to the solutions of the Boltzmann equations]^ 



3 As discussed in Section [3J we use an approximate solution of the Friedmann equation, H — a/a, with 



the scale factor a(t) given in Eq. (36). 
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Figure 5: Contour plots of the reheating temperature Trm as a function of the parameters fhi and Mi 
for the three different choices of vb-l listed in Eq. ( 102 1. The outcome of the Boltzmann equations Trh, 
calculated according to Eq. (105), is compared with three different estimates T^jj, T^jj and T^jj which 



are respectively defined in Eqs. (106), (C.8) and (C.9) . The contour labels as well as the background 



colours indicate the numerical values of Trh ■ Going to smaller values of mi the T RH - and T^h contours 

(2) 

approach the corresponding Trh contours from below. With respect to the Trh contours the T RH - 
contours are shifted upwards by approximately A log 10 Trjj 



0.04, cf. Eq. (C.ll). 
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For nonrelativistic Majorana neutrinos Nf, instantaneous energy transfer and H being 
determined from the Friedmann equation, the reheating temperature is given by 



In our scenario we do not meet either of these conditions which is why the simple estimate 
T RH has to be augmented with several corrections in order to properly reproduce the 
outcome of the Boltzmann equations: 

• Being produced in S decays, the nonthermal neutrinos all carry initial energy 
ms/2 — 150Mx- For most of the time they are, hence, highly relativistic such 



that their decays occur at an effective rate (cf. Eq. (105)). Replacing T^, by 



Tfj results in an approximation T RH for the reheating temperature; 

At the time when Trh is evaluated a large fraction of the energy density still resides 
in nonthermal neutrinos. Taking into account that only a part of the total energy 
density at t = tun contributes to Trh yields an approximation T RH ; 

The fact that our approximation for the Hubble rate H = a/a does not fulfill the 
Friedmann equation exactly introduces a final imprecision which effectuates the 

(2) 

remaining small deviation of T RH from the actual reheating temperature Trh- 



We refer the interested reader to Appendix O where the reconstruction of the numerical 
result Trh starting from the simplest estimate T RH is discussed in greater detail. 

The reheating temperature Trh obtained from the Boltzmann equations for the three 



values of v b-l m Eq. ( 102 ) is presented in Fig.|5jtogether with the different approximations 
T RH as a function of the neutrino mass parameters rh\ and M\. Notice that the behaviour 



of T RH is determined by the width ex mi M\ (cf. Eq. (12)), which is independent of 



vb-l- The correction due to the time dilatation factor mainly depends on the ratio of the 



Majorana neutrino decay width Eq. (12) and the S boson decay width Eq. (22) 



T° Nl m lV 2 B _ L 

~W K • ( 107 ) 

For ^> Tg, the bulk of the nonthermal neutrinos decaying at t = tRH is produced 
just shortly before and is therefore relativistic. On the other hand, for <C T° s , most 
of the nonthermal neutrinos decaying at the reheating time are nonrelativistic. For fixed 
vb-Li this correction turns out to be marginal for the smallest effective neutrino masses 
rh\ and the largest Majorana neutrino masses M\. The correction increases with the ratio 
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in Eq. (107) becoming larger. Its maximum is given by the flavour model, ^/iMijms 



^ l /lJl5Q. The related correction corresponding to the overestimation of the energy density 
of radiation at t — t R n has the same dependence on parameters. Finally, the mismatch 
between the Hubble rate and the exact solution of the Friedmann equation only slightly 
modifies the reheating temperature. All in all, the global effect of these corrections is to 
increase (decrease) the dependence of the reheating temperature on Mi (mi). 

In each of the three panels of Fig.[5j corresponding to the three different values of 
vb-l, the values of rhi and M\ respectively span four orders of magnitude allowing for 
reheating temperatures ranging over five orders of magnitude. Reheating temperatures 
as small as T RH ~ 10 5 GeV are obtained for the lowest decay rates in association with the 
smallest initial false vacuum energy density, i.e. for the minimal values of vb-l, ™i and 
Mi. Conversely, reheating temperatures as large as Trh — 10 12 GeV are obtained for the 
maximal values of vb-l, fhi and Mi. 

Lastly, we observe that the region where the reheating temperature exceeds the Ma- 
jorana neutrino mass significantly shrinks when going from the simplest approximation 
Tftjj to the results of the Boltzmann equations T RH . As for the former, T R ^ > Mi for 
rhi ^2x 10 _3 eV, independent of Mi, while in the latter case Trh > Mi is only accom- 
plished for the largest values of rhi and Mi. The reasons for this relative decrease in the 



reheating temperature were already mentioned following Eq. (106): the longer neutrino 



lifetimes due to their relativistic nature and the overestimation of radiation energy in de- 
riving T r °jj. As the strength of the washout processes during the generation of the baryon 
asymmetry crucially depends on the ratio of temperature T and neutrino mass Mi we 
expect the efficiency of leptogenesis to severely drop in the region T RH > Mi . 

5.2 Baryon asymmetry 

The baryon asymmetry that is generated for a given choice of input parameters follows 



from the respective solutions of the Boltzmann equations according to Eq. (91). In this 
section we shall discuss in turn the contributions it receives from the decays of the non- 
thermal and thermal neutrinos. Our main results are displayed in Fig.[6] which presents 



the baryon asymmetry for the three values of the B — L breaking scale (cf. Eq. (102)) as 
function of the neutrino mass parameters rhi and Mi. In each panel of Fig.[6| we indicate 
the regions in which leptogenesis from the decay of nonthermal (light green) and ther- 
mal neutrinos (gray-green) successfully produces the observed baryon asymmetry. Notice 
that in some regions of parameter space (dark green) both variants of leptogenesis man- 
age to overcome the observational bound individually while in others (white) r]° B hs only is 
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exceeded after taking the sum of the two contributions. 

The blue solid lines in Fig. [6] separate the parameter regions in which leptogenesis is 
respectively driven either by the decay of the nonthermal or the thermal neutrinos. In the 
viable regions of parameter space the nonthermal contribution to the baryon asymmetry 
typically represents the clearly dominant one. As we have checked numerically, the gen- 
erated asymmetry in that case can be reconstructed to good approximation by assuming 
that at t — t\ the energy density of the nonthermal neutrinos is almost instantaneously 
converted into radiation. Nonthermal neutrinos of average energy ef^ (cf. Eq. (|84[)), that 
rapidly decay around £ = t±, lead to an asymmetry (cf. [37] ) 



q rp 

rapid 1-7 u ± 



(10* 

t=ii 



-JVi 

We emphasize that neglecting the relativistic motion of the nonthermal neutrinos, i.e. 



employing simply the mass M\ instead of the full energy per particle e s N in Eq. (108) 



would entail an asymmetry proportional to the temperature, rf^ lA oc T(ti). Such an 
estimate fails to reproduce our results except for some accidental points in parameter 
space Q 

For fixed Vb-l and rh\ < (9(10~ 3 )eV the nonthermal baryon asymmetry does not 
depend on rh\ anymore. This observation can be easily understood in terms of the Boltz- 



mann equation (80a) for N§_ L . For very small effective neutrino masses, the washout 
processes become inefficient, leaving us only with the production term. The size of the 
final asymmetry then only depends on the maximal abundance that can be reached 
in the course of S boson decay, which is reminiscent of standard thermal leptogenesis in 
the weak washout regime. Since the collision operator accounting for the production of 
nonthermal neutrinos through the decays of S bosons as well as the CP parameter e\ are 
solely controlled by the Majorana neutrino mass, the resulting baryon asymmetry ends 
up being exclusively determined by M%. 

Increasing Vb-l for fixed neutrino masses rh\ and M\ reduces the produced baryon 



asymmetry. This is due to several effects whose influence is apparent in Eq. ( 108 ): On the 



one hand a higher B — L breaking scale implies a larger relativistic correction resulting 
in a smaller effective decay rate rjy , on the other hand it leads to a faster Hubble 
expansion. The former increases efj ~ (Mi/ En^ 1 Mi and delays the neutrino decays 



such that Eq. (108) needs to be evaluated at a later time t\ corresponding to a smaller 



temperature T. The faster Hubble rate H reinforces the drop-off in the temperature. We 

7 This actually happens in Ref. [3] in which 77^ pld is calculated for the corresponding set of parameter 
values given in Tab. [2j 
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Figure 6: Contour plots of the baryon-to-photon ratio t\b 'as defined in Eq. (91) as a function of the 



parameters rh\ and Mi for the three different choices ofvs-h listed in Eq. ( 102 1. In the light green (gray- 
green) region leptogenesis through the decay of nonthermal (thermal) Majorana neutrinos successfully 
reproduces the observed baryon asymmetry rfg s . In the red region the total asymmetry is not able to 
overcome the observational bound. The blue line separates the regions in which each leptogenesis variant 
is the dominant one. 
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may reformulate this argument in terms of the rates Tf^ and H by saying that a smaller 
ratio Fff /H reflects a lower efficiency of the nonthermal neutrino decays. From this point 
of view, the generation of the asymmetry struggles to keep pace with the expansion of 
the universe resulting in a more diluted asymmetry. 

At values of rh\ larger than O(10~ 3 ) eV resonant £H -H- £H scatterings that wash out 
the generated asymmetry at a rate (cf. Eq. (94)) dramatically decrease the efficiency 
of nonthermal leptogenesis: For small rh\ we have Trh ^ M (cf. Fig. [5]) and the 
production of on-shell N\ neutrinos out of the thermal bath is Boltzmann suppressed, 

T«M 1 : T W = Jlr^ oc e-W Y% « T° Ni . (109) 

But as rh\ becomes larger, the reheating temperature approaches M\ and the final asym- 
metry is depreciated due to washout. For given rh\ and M\, increasing vb-l results in 



a decrease of the reheating temperature (cf. Section 5.1). In particular, this reduces 
the region in parameter space where Trh > Mi, consequently extending the region in 
which nonthermal leptogenesis can successfully proceed without being much affected by 
washout. 

The decay of the nonthermal neutrinos is not the only mechanism by means of which 
the baryon asymmetry is generated in our scenario. It also receives a contribution rf^ from 
the decays of the thermally produced neutrinos N± . If M\ is sufficiently large, rjg can 
exceed the observed baryon asymmetry on its own. Its behaviour in parameter space is 
similar to the one of standard thermal leptogenesis with vanishing initial neutrino abun- 
dance. However, it is important to note that standard thermal leptogenesis differs from 
our thermal mechanism in the sense that in the former the thermal bath out of which the 
Majorana neutrinos are produced is assumed to have an independent origin (e.g., inflaton 
decay) and the initial temperature usually is taken to be arbitrarily high. By contrast, the 
generation of our thermal asymmetry rjg is tightly coupled to the dynamics of reheating 
in the course of the nonthermal neutrino decays. In standard thermal leptogenesis the 



CP asymmetry e\ (cf. Eq. (14)) as well as the evolution of the N± and B — L abundances 
are controlled by the neutrino mass parameters fhi and M\. To guarantee successful lep- 
togenesis, Mi is constrained to be at least (9(10 9 ) GeV if rh\ is fixed at rh\ ~ 10 -3 GeV. 
Effective neutrino masses rh\ different from that result in larger bounds on Mi. 

In our scenario, the values of Mi above which thermal leptogenesis is efficient are 
comparatively one to two orders of magnitude larger. Tab. [3] summarizes the respective 
bounds on Mi for the three different B — L breaking scales together with the corresponding 
values of rh\ and Tr#. The fact that now Mi has to be much larger than (9(10 9 ) GeV 
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Panel Vb—l [GeV] fhi [eV] M x [GeV] TW [GeV] T RH /M 1 



1 


3.4 x 10 12 


3.2 x 10" 3 


1.0 x 10 10 


6.4 x 10 9 


0.63 


2 


5.8 x 10 13 


5.2 x 1(T 3 


2.1 x 10 10 


7.9 x 10 9 


0.38 


3 


1.0 x 10 15 


1.6 x 1(T 2 


8.5 x 10 10 


2.0 x 10 10 


0.23 



Table 3: Parameter points in the three panels of Fig.\Qcorresponding to the lowest possible values of Mi 
for which the decay of the thermal neutrinos suffices to reproduce the observed baryon asymmetry. The 
values for the reheating temperature follow from Fig. [5l 



finds its origin in the interplay between the specific reheating process at work and the 
temperature dependence of thermal leptogenesis: First of all, in the discussion of the 
decay of the nonthermal neutrinos we saw that the temperature is bounded from above 
to prevent complete washout of the asymmetry. The same holds for thermal leptogenesis; 
but in this case the temperature also must not be too low in order to ensure an efficient 
neutrino production from the thermal bath. Consequently, as a compromise between 
very small (T <C Mi) and very large (T ^> Mi) temperatures, thermal leptogenesis is 
most efficient at T ~ M\ (cf. Tab. |3])j^] Second, as for standard thermal leptogenesis, 
our thermal mechanism prefers an intermediate value of fhi. Taking fhi to large values 
increases the strength of the washout processes. Small fhi results in a low temperature and 
a small neutrino decay rate T° Ni such that the Nf production becomes suppressed. When 
asking for the lower bounds on Mi we thus have to look for the smallest values of Mi for 
which the condition T ~ Mi holds and fhi is neither too large nor too small. In contrast 
to standard thermal leptogenesis, in our scenario the accessible temperatures now also 
depend on the mass parameter Mi. As can be seen from Fig.[5j the considered reheating 
process simply does not manage to satisfy the condition T ~ Mi for Mi ~ 10 9 GeV 
without entering the strong washout regime. Instead, Mi has at least to be as large as 
indicated in Tab. [3] to avoid too large values of fhi while still fulfilling T ~ Mi. 

Comparing the three points in Tab. [3] we note that the ratio of the reheating temper- 
ature Tru to Mi decreases as vb-l becomes larger. The production of thermal neutrinos 
is, consequently, less efficient for high vb-l- This is, however, compensated for by the 
increase in the CP asymmetry parameter ei for heavier N\ neutrinos (cf. Eq. (15)). Like- 
wise, the corresponding effective neutrino masses fhi increase when going to larger B — L 



8 Note that our scenario also differs from standard thermal leptogenesis because we only consider 
decays and inverse decays. We have checked that including AL = 1 and AL = 2 scatterings would 
enforce the production of Majorana neutrinos for mi <3x 10~ 3 eV as it is the case in standard thermal 
leptogenesis |33j . resulting in a slight expansion of the allowed region in the weak washout regime. 
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breaking scales. This effect is based on the fact that for fixed mi and Mi an increase in 
Vb-l entails a drop in the temperature. The factor representing the Boltzmann suppres- 



sion in the washout rate IV (cf. Eq. (109)) then becomes smaller which enables one to 



raise the neutrino decay width T° Ni by increasing m\. 

Standard thermal leptogenesis predicts a final baryon asymmetry of 

Vb = ^-c sph eiK f (mi) . (110) 

where the final efficiency factor Kf only depends on rh\. For rh 1 > 10 _3 eV it may be 
parametrized as [33] 

o /10- 2 eV\ u , N 

K f (mi = 2 x 10~ 2 — . Ill 

V mi J 

Combining Eqs. (110) and (111) with Eq. (15), one finds that rf^ evolves as rf^ oc Mi. 



This is exactly the behaviour of the total baryon asymmetry one observes in the regions 
where the thermal contribution dominates over the nonthermal one, i. e. the regions on the 
right-hand side of the blue lines in Fig.[6j As the number density of nonthermal neutrinos 
usually exceeds the number density of thermal neutrinos at the time the asymmetry is 
created, the relative size of the two asymmetries 7/J> and rf^ is controlled by the efficiency 
of the nonthermal mechanism. Only when the nonthermal asymmetry is suppressed due 
to efficient washout, the baryon asymmetry due to the decay of the thermal neutrinos has 
a chance to dominate. 

In conclusion, it is remarkable that leptogenesis through the decay of the nonthermal 
Majorana neutrinos is able to widely extend the region in parameter space in which the 
observed baryon asymmetry can successfully be reproduced. For the lowest B—L breaking 
scale vb-l = 3.4 x 10 12 GeV, Majorana neutrinos as light as Mi ~ 10 8 GeV are sufficient 
to generate the observed baryon asymmetry. 

5.3 Gravitino dark matter 

Having discussed leptogenesis on its own in the last section, we now ask for the regions 



in parameter space where both conditions of Eq. (101) are satisfied, i.e. in which we 



obtain gravitino dark matter along with a sufficient baryon asymmetry. As outlined in 



Section 3.2.6 the thermal production of gravitinos is controlled by three parameters: the 
gravitino and gluino masses mg and rrig as well as the temperature T. The latter is 
determined by the reheating process, T RH = T rh (vb-l, m i, M), such that Q^h 2 , the 
present contribution from gravitinos to the energy density of the universe, depends on 
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all free parameters of our scenario. For each point in parameter space the respective 



solutions of the Boltzmann equations allow us to calculate Qgh 2 according to Eq. (98). 



By imposing the condition that gravitinos be the constituents of dark matter we can then 
derive relations between neutrino and superparticle masses. For instance, if we fix the 
gluino mass at 800 GeV, 



n d h 2 (v B . L , Aft, mi, mg, mg)\ m _ =moGey 
we can solve for Mi as a function of vb-l, m\ and mg, 

Mi = Mi(w B _ L , mi,mg) . 



(112) 



(113) 



We consider those choices of the parameters vb-l, m i an d mg as viable, which actualize 
gravitino dark matter for Mi values that are accessible in the context of the flavour model 



(cf. Eq. (7a)) 



Mi(v B -l, mi, mg) < rfvB-L ■ 



(114) 



Furthermore, applying Eq. (113) to the results of Section 5.1 allows one to trade the Mi 
dependence of the reheating temperature Trh for a dependence on mg, 



T RH = T rh (v B -l, Mi, mi) ->■ T rh (v B -l, m h mg) . 



(115) 



Similarly, Eq. (113) can be used to translate the bounds on mi and Mi shown in Fig. [6 
that were obtained by requiring successful leptogenesis into bounds on m x and mg, 



i] B (v B -L,M 1 ,m 1 ) ->■ i] B (v B ~L,m 1 ,m 5 ) >r]° B bs . 



(116) 



The parameter points we are after, i.e. the points at which the baryon asymmetry is 
accounted for by leptogenesis and gravitinos constitute the dark matter, now correspond 



to those values of vb-l, mi and mg that satisfy the two conditions in Eqs. (114) and 



(116) simultaneously. On the basis of our numerical study of the Boltzmann equations 



we are able to identify the regions of interest in parameter space: Fig. [7] presents our 
results in combination with the associated values of Mi (cf. Eq. ( 113| )), Fig. [8] features 



the related reheating temperatures (cf. Eq. (115)) instead. Again, both figures consists 
of three panels each that respectively take care of the the three different B — L breaking 



scales specified in Eq. (102). We refer the interested reader to Appendix D which gives a 



detailed account of how Figs. [7] and [8] can be reconstructed by means of simple analytic 
expressions and with the aid of our numerical findings for Trh and r\B- 
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Figure 7: Gravitino mass range consistent with gravitino dark matter (cf. Eq. (112)) and successful 
leptogenesis (cf. Eq. ( |116[ )j depending on the effective light neutrino mass. The contour lines refer to 



the neutrino mass Mi (cf. Eq. (114)) as a function of rh\ and mg such that the gravitino abundance is 
ilgh 2 = 0.11. In addition to the colour code introduced in Fig.[fi[ the black region in the upper-left panel 
represents the Mi values that are not allowed by the flavour model. 
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Notice that we also consider gravitino masses almost as large as the gluino mass, 
m c — ^00 GeV while m g = 800 GeV. Imposing gaugino mass unification at the GUT scale 
would, however, forbid such a nearly degenerate superparticle spectrum. The running of 
the renormalization group equations would then imply a gaugino mass relation M3/M1 ~ 
5.9 at low energies. Given that it is the lightest supersymmetric particle, the gravitino 
would have to be lighter than the bino resulting in an upper mass bound of mg < 140 GeV. 

If we were to select a gluino mass other than m g = 800 GeV all values of in Figs. [7] 
and [8] would have to be rescaled while vb-l, Mi and rh\ could remain unchanged. This 
follows from the fact that the gravitino abundance in Eq. (112) can be kept constant by 
compensating a change m g — > am g in the gluino mass by a change mg — » 6mg in the 
gravitino mass without altering the reheating temperature Trh- As long as mg is much 
smaller than m§, the factor b simply corresponds to a 2 . The general relation between a 
and b is discussed in Appendix [Dj To sum up, thanks to this relation the results presented 
in Figs. [7] and [8] can be generalized to different gluino masses by correspondingly relabeling 
the gravitino axis. 

As a general trend in Figs. [7] and [8] we observe that for fixed rh\ and mg < 230 GeV 
the neutrino mass M\ and the reheating temperature Trh continuously become larger 
when increasing the gravitino mass. For > 230 GeV this behaviour is reversed: The 
M% contours in Fig. [7] bend over as soon as ~ 230 GeV is exceeded. In Fig. [8] the 
temperature remains rather constant at Trh ~ 5 x 10 9 GeV for 150 GeV < < 400 GeV. 
Beyond > 400 GeV it begins to decrease againj^] The physical origin of these two 
regimes can be traced back to the rate = 7<g/ng (cf. Eqs. (85) and (86)) at which 
gravitinos are created from the thermal bath (cf. Eq. ( D.6[ )) [1~4"][T5] . 

r g = r g (T,m g ,m,)(x ^l + ^|Pj . (117) 



In the regime mg <C m g (T) the second term in Eq. (117) is the dominant one and it is 
mainly the goldstino part of the gravitino, i.e. its components with helicity ±|, that is 
produced. A larger gravitino mass then implies a smaller rate necessitating a stronger 
reheating in order to still generate the right abundance. Correspondingly, the neutrino 
mass Mi also has to increase to bring about the higher temperature. Evolving a gluino 
mass of 800 GeV from the electroweak scale to a temperature T ~ 5 x 10 9 GeV results 
in a high-scale mass of m g (T) ~ 400 GeV. Because of that, is dominated by the first 
term in Eq. (117) from mg ~ 400/v^3GeV ~ 230 GeV onwards. This means that, for 



9 Cf. the contour corresponding to Trh — 10 9 ' 5 GeV reentering the second panel of Fig. [8] at mi 
lCT 5 eV and m d ~ 700 GeV. 
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such large values of mg, primarily the transverse degrees of freedom of the gravitino, 
i.e. its components with helicity ±|, are excited. The production rate then becomes 
independent of mg turning into a function of the temperature T only. In such a case the 
final gravitino abundance f2g/i 2 simply scales linearly with mg (cf. Eq. (D.l)). Hence, 
larger gravitino masses have to be balanced by smaller reheating temperatures to keep 
Q(jh 2 fixed. This explains the decrease in Trh and Mi at very large gravitino masses. 

On the other hand, varying rhi at constant mg has almost no effect on the reheating 
temperature, which is expected since Tg inherently is a function of T, mg and trig. As 
each gravitino mass is associated with an appropriate rate Tg, the choice of mg already 
implies a unique reheating temperature T RH f» T RH (m^), independent of the underlying 
neutrino parameters (cf. Fig. 12 in Appendix D). Meanwhile, the neutrino mass Mi 
becomes smaller when increasing rh\ in order to ensure that T R h remains approximately 
constant for fixed mg (cf. Fig. [5]). 

In Figs.[7]and[8]we also indicate the regions in parameter space that are not compatible 
with our scenario because either of the two conditions in Eqs. (114) and (116) is not 
satisfied. Bounds coming from the flavour model (cf. Eq. (114)) only show up for vb-l — 



3.4 x 10 12 GeV: The requirement that M 1 be smaller than 1.1 x 10 10 GeV (cf. Eq. gggp ) 
excludes gravitino masses larger than 35 GeV for rhi = 10~ 5 eV. At rhi = 10~ 3 eV it rules 
out masses in the range between 120 GeV and 430 GeV and for rh 1 > 1.8 x 10 _3 eV it 
does not constrain mg any longer at all. In the case of the two other choices for vb-l the 
respective flavour bounds on M\ are never reached because the corresponding reheating 
temperatures are too high. Demanding a sufficient baryon asymmetry (cf. Eq. (116)) 
yields lower bounds on mg in the weak washout regime and limits the maximal value of 
rh\. Notice that these bounds are in one-to-one correspondence with the constraints on M\ 
and rhi i n Fig.[6j For instance, at small rhi the gravitino mass can only decrease as long 
as Mi is large enough so that the observed baryon asymmetry is reproduced. Similarly, 
at large rhi the sharp drop-off in the efficiency of leptogenesis due to stronger washout 
limits the viable range of rh\. In Tab. [4] we present the smallest gravitino masses that are 
accessible for certain representative values of rh\. As in Fig. [7]the contour lines of constant 
Mi slightly fall off with decreasing rhi, we find the lowest bounds on mg at rhi = 10 -5 eV. 
For weak washout higher B — L breaking scales lead to tighter bounds on mg, just as it is 
the case for the neutrino mass Mi (cf. Fig. [6]). In the strong washout regime we encounter 
the opposite behavior. Here, the contour line corresponding to i]b = rj B hs , which separates 
the allowed and excluded regions in parameter space, steeply rises. In Section |5.1| we 



argued that the larger the value of vb-l the later this rise sets in when increasing mi 
(cf. Eq. (109)). Therefore, the tightest bounds on mg are now obtained for low B — L 
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Panel v B -l [GeV] / m x [eV] 10" 5 10" 4 lO" 3 10~ 2 10' 



1 3.4 x 10 12 


180 MeV 


360 MeV 


870 MeV 




2 5.8 x 10 13 


570 MeV 


1.2 GeV 


2.5 GeV 


70 GeV 


3 1.0 x 10 15 


1.9 GeV 


3.7 GeV 


6.2 GeV 


9.2 GeV 



Table 4: Lower bounds on the gravitino mass according to Figs, for the three different choices 



of vb-l listed in Eq. ( 102 ) and five different values of rk\. A dash ( — ) indicates that leptogenesis is not 
efficient enough to produce the observed baryon asymmetry as long as the requirement of gravitino dark 
matter is kept. 



breaking scales. 

We also note that for rh\ = 0.1 eV it is not possible to produce the observed baryon 
asymmetry while sticking to the assumption of gravitino dark matter, independent of the 
value chosen for Vb-l- As rh\ is bounded from below by mi, the smallest eigenvalue 
of the standard model neutrino mass matrix, this observation opens up the possibility 
of falsifying our proposed scenario in future neutrino experiments. The measurement 
of a light neutrino mass of 0.1 eV, combined with the known differences of the light 
neutrino masses squared, would imply that rh\ > 0.1 eV, thereby ruling out our mechanism 
of entropy production. Likewise, any lower limit on the absolute neutrino mass scale 
coming from, e.g. cosmological observations would restrict the allowed range for the 
gravitino mass. A determination of the gravitino mass on the basis of cosmic gamma-ray 
observations or decays of the next-to-lightest-superparticle (NLSP) in collider experiments 
could, in turn, constrain the neutrino mass spectrum. 

In standard thermal leptogenesis the reheating temperature has to be at least T RH > 
10 9 GeV, independent of the initial conditions, to guarantee a successful generation of the 
baryon asymmetry [33J. Together with the lower bound on the gluino mass imposed by 
collider searches, this constrains the gravitino mass to lie in the range mg ~ 10 100 GeV 
in order to be compatible with the observed dark matter abundance. By contrast, the 
present scenario allows for a much broader range of gravitino masses since the reheating 
temperature can be significantly lower than in the case of thermal leptogenesis. As appar- 
ent in Fig.[8j Trh can decrease down to values of O(10 7 ) GeV if fh\ and M\ are chosen such 
that the nonthermal neutrinos decay extremely slowly. We note that for such reheating 
temperature production of gravitinos from inflaton decay is usually negligible [3H]. This 
paves the way for gravitino masses as small as 200 MeV (cf. Tab. [4]). On the other hand, 
as our scenario can as well accomodate neutrino masses M\ of order 10 10 -j- 10 11 GeV, 
the gravitino can also be almost as heavy as the gluino, mg ~ few x 100 GeV. Note, 
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however, that requirements such as gaugino mass unification will further constrain the 
superparticle spectrum. 

Finally, we would like to point out that low gravitino masses have interesting con- 
sequences for the decay of the next-to-lightest superparticle (NLSP). If -R-parity is con- 
served, the lower bound on the mass of stable gravitinos from standard thermal leptoge- 
nesis, mg > 10 GeV, implies a long NLSP lifetime which could jeopardize the success of 
primordial nucleosynthesis (BBN). 

A study of general neutralino NLSPs with gravitino LSP has been performed in [39] . 
In this work, lower mass bounds for different NLSP types have been extracted from the 
bounds on the decay of heavy neutral particles during BBN for given gravitino masses. 
Assuming a gravitino of 10 GeV, it turns out that the predictions of primordial nucle- 
osynthesis are not affected for pure bino, pure wino and mixed gaugino-Higgsino neu- 
tralino NLSPs that are heavier than 3, 0.8 and 1 TeV, respectively. These bounds are 
now significantly softened for the light gravitino masses which can occur in our sce- 
nario: For mg = 200 MeV, pure bino, pure wino and mixed gaugino-Higgsino neutralino 
NLSPs as light as 800, 150 and 200 GeV, respectively, are in agreement with the BBN 
bounds. These lower masses are interesting for two reasons: First, they allow a smaller 
scale of supersymmetry breaking and second, they can be probed more easily at the 
LHC. For a charged NLSP like a scalar r-lepton, its lifetime has to be sufficiently short, 
tnslp < 10 3 -10 4 s [IKE], which typically requires m d < 1 GeV for m NLSP = 0(100) GeV. 

Remarkably, the small gravitino masses required from such constraints on NLSP decays 
can be accommodated in our scenario. We thus emphasize that reheating through the 
decays of heavy neutrinos provides a way to alleviate the existing tension between the 
generation of the baryon asymmetry, stable gravitino dark matter and BBN constraints. 

6 Summary and outlook 

We have studied the production of the entropy of the hot early universe through the 
decays of heavy Majorana neutrinos. As an example, we considered the decay of a false 
vacuum with unbroken B — L symmetry. Tachyonic preheating leads to a state whose 
energy density is dominated by nonrelativistic S bosons, the Higgs boson associated with 
spontaneous B — L breaking, with a subdominant admixture of radiation. Subsequent pro- 
duction of the lightest heavy Majorana neutrinos N\ from S decays and from the thermal 
bath leads to a phase whose energy density is dominated by N\ neutrinos. Their decay 
then produces all entropy of the hot thermal universe, together with the baryon asym- 
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metry via a mixture of nonthermal and thermal leptogenesis. Simultaneously, thermal 
production of gravitinos accounts for the observed dark matter. 

We have studied the time evolution of this system by means of a set of Boltzmann 
equations for distribution functions which take into account the differences between ther- 
mally and nonthermally produced Ni neutrinos. Details of the initial state are rather 
unimportant for the final result. What matters is the existence of a phase dominated by 
Ni neutrinos. Their lifetime determines the reheating temperature that is crucial for the 
final gravitino abundance and hence for the produced amount of dark matter. 

Our analysis has been based on a flavour model that allows to vary the key parameters 
for leptogenesis, the neutrino masses Mi and rhi, over a wide range, consistent with the 
masses and flavour mixings of charged leptons and neutrinos. The final baryon asymmetry 
and the dark matter abundance have been calculated in terms of several parameters of 
the Lagrangian, independent of initial conditions: the scale vb-l of B — L breaking, the 
heavy Majorana neutrino mass M l5 the effective light neutrino mass fhi, the gravitino 
mass itlq and the gluino mass rrig. For generalisations of the flavour model the produced 
matter-antimatter asymmetry and dark matter can depend on further parameters. 

Particularly interesting is the resulting relation between the lightest neutrino mass 
and the gravitino mass. For instance, for a typical gluino mass of 800 GeV and a light 
neutrino mass of 10~ 5 eV the gravitino mass can be as small as 200 MeV, whereas a lower 
neutrino mass bound of 0.01 eV implies a lower bound of 9 GeV on the gravitino mass. 
The measurement of a light neutrino mass of 0.1 eV would falsify the proposed mechanism 
of entropy production. These results provide an important connection between collider 
searches for superparticles and neutrino mass determinations in laboratory experiments 
and by cosmological observations. Measurements of the absolute neutrino mass scale and 
superparticle masses consistent with our predictions would provide important indirect ev- 
idence for the origin of entropy, matter and dark matter and for the maximal temperature 
of the hot thermal universe. 

We are planning to extend our analysis in several directions: Boltzmann equations 
for the superpartners of S bosons and heavy Majorana neutrinos have to be incorporated 
in a fully consistent calculation of baryon asymmetry and gravitino dark matter. Further 
important questions concern the connection with inflation and the possible relation 
between B — L breaking and supersymmetry breaking. 
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A Conventions for the Boltzmann equations 

The Boltzmann equation for a particle species ip describes the time evolution of the 
distribution function /^(t,^^,^) in the one-particle phase space [37], where is 
defined such that f^d^^ gives the average number of ip particles in the phase space 
volume G?$^ at time t. Imposing homogeneity and isotropy of spacetime, only depends 
on the absolute value p of the 3-momentum p as well as on time t. The Boltzmann 
equation for ip particles in the Friedmann-Lemaitre framework then reads 

CU(t,p)=(J- t -Hp^jU(t,p)=C^, (A.l) 

where with £ we denote the Liouville operator, and with the collision operator. The 
latter keeps track of changes in due to inelastic interactions, and may be decomposed 
into contributions from decays and scatterings, 

ij... a ij... 

The operators are obtained from quantum mechanical transition probabilities in- 
tegrated over the multi-particle phase space 

C^ab...^i 3 ...) = -±— [ dn^|a,6,...^,j,...)(27r) 4 5( 4 )(EPo Ut -EPin) (A.3) 

x /,/•,.. (1 ± U) (1 ± /„ ) (1 ± f b ) ... M {ij... > cab...) 2 
- /,,/•„/,,.. (1 ± fi) (1 ± ./• ) ... M [cab... - •> ij...) 2 } , 

where is the number of internal degrees of freedom of ip and dll subsumes all Lorentz 
invariant momentum space elements dp = (2ir)~ 3 d 3 p/2E along with a statistical factor S 
that prevents double counting in the case of identical particles 

dll(il)\a,b, ...) = S(ip,a,b, ...)dp a dp b ...dp i dp j ... . (A.4) 

The amplitudes squared |-M| 2 are understood to be summed over all internal degrees of 
freedom. Since they are expected to yield only small corrections [34], the Bose enhance- 
ment (1 + /) and Pauli blocking (1 — /) quantum statistical factors related to boson and 
fermion production respectively, are neglected in this work. Their influence may partly 
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be canceled by other quantum corrections, like off-shell effects [12]. Additionally, the 
operators can be split into direct and inverse processes 



C^{ipab... -H- ij...) = C^iij... — > ipab...) — C^ijpab... — > ij...) . 



(A.5) 



If the ip particles are in kinetic equilibrium, the integration of Eq. (A.l) over the ip 
phase space leads to a Boltzmann equation for the ip number density 



(A.6) 



y- 



where n^, and the interaction densities 7 are defined as 



j(i[)ab... <-> ij...) 



(2tt)' 



d 3 p C^(ipab... ^ ij...) . 



(A.7) 
(A.8) 



The Boltzmann equation Eq. (A.6) can alternatively be written as an equation for the 



comoving number density A^, = a 3 n^ as a function of the scale factor a 



da 



aH—N^ = a 3 



(A.9) 



B Phase space distribution of thermal neutrinos 



When deriving the Boltzmann equation (69) for thermal neutrinos in Section 3.2.3, we 



asserted that these are approximately in kinetic equilibrium (cf. Eq. (68)). Given a 
numerical solution for T(a) for a specific choice of input parameters Vb-l, fn>i and M 1; we 
can check the self-consistency of this approach by comparing our approximate distribution 



function proportional to fff with the exact expression in Eq. (67). In this appendix we 
perform such a comparison for the exemplary parameter point discussed in Section [4] (cf. 
Tab. [2]). To begin with, we introduce the following ratio 

fl(t, P ) 

yjg(*>p)' 



R T Nl (t,p) 



t 




J dt' exp 


-Mxl^y* dt"E^(f)j 


t 2 





Mi r0 m?, P ) 

J- T 



E Nl (t')^^f%(t,p) 



(B.l) 
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Figure 9: Ratio i?^ = fjf 1 /fj? normalized to its mean value u.r as a function of the neutrino momentum 
p at different times, cf. Eqs. ( |B.2[ ) and (B.3) as well as Tab. [5] The input parameters Vb-l, Si and 
Afi are chosen as in Section [I] (cf. Tab. [2|). The respective values of the scale factor are discussed in 



Section 4.2 The white and gray bands separate the different orders of magnitude. 



with /jy taken from Eq. (67), and determine its momentum dependence at different times. 
A momentum independence of Rjj would then reflect an exact kinetic equilibrium. In 
this case the equation = R% can easily be integrated over phase space yielding 



R T Nl = R T Nl (t) O (t) = R T Ni {t)N% (t) & f T Nl (t, p) = (t, P) ■ (B.2) 



A convenient measure for the deviation from kinetic equilibrium at a given time is the 
coefficient of variation cr = <Tr/ {J,r, i.e. the standard deviation (Xr of Rj? in relation to 
its mean value /ir with respect to an appropriate momentum interval Ap, 

Mt) = (R T Nl ) P > = ((« 1 ) 2 ) p - « 1 >') 1/2 ( B - 3 ) 

where (•) is defined as 

Ap 

1 



{f) ? {t) = Ap I d PM>P)> (R4) 



and Ap, by convention, is always chosen as 



/£(t,Ap) = l(r 4 x/*[(t,0) , (B.5) 
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such that the relevant range of momenta is covered. 

We compute R% , fiR, ctr, cr and Ap for six representative values of the scale factor, 
and summarize the corresponding results in Fig. [9] and Table [5j At early times small 
momenta are much more frequent than in kinetic equilibrium and states with large mo- 
menta are underpopulated. For extreme momenta, Rjf can become ten times as large or 
small as its mean value /j,r. As time goes on, this tilt in R% , however, disappears and 
for a ~ 12000, kinetic equilibrium is eventually reached. On average Rj^ deviates from 
Hr not more than one order of magnitude and, from this perspective, the approximation 
of kinetic equilibrium may be regarded as justified. The steady convergence to kinetic 
equilibrium is also reflected in the behaviour of the coefficient of variation cr which starts 
out at a value of cr ~ 0{1) and decreases to cr ~ (9(10~ 4 ). For other choices of the 
model parameters we expect the Nf phase space distribution to behave similarly. 

Scatterings of the thermal neutrinos involving standard model quark pairs such as 
Ni£ -H- qu, Niu -H- £q and N±q <-> £u, speed up the equilibration of the neutrino dis- 
tribution function [33]. This results in a larger abundance of thermal neutrinos at high 
temperatures. On the other hand, scatterings also tend to increase the efficiency of 
washout processes such that, after all, their impact on the generated thermal asymmetry 
is negligible for our purposes. 

There are two main reasons why Rjj is not flat from the beginning: The first is 
directly related to the momentum dependence of the production and decay terms in the 



Boltzmann equation (69) for thermal neutrinos. In both terms the effective rate T Ni 
comes weighted with the inverse time dilatation factor (Mi/E^). It is, thus, larger 
at smaller momenta which is why initially, when <C f^, neutrinos with smaller 
momenta are overproduced. Once has overshot the decay term dominates, again 
preferably changing the abundance of low-momentum states. This interplay between 
production and decay is balanced such that R% is eventually flattened out. A numerical 



investigation of the different factors in the integrand of Eq. (B.l) confirms this simple 
argument. The second reason is the intermediate stage of reheating between the phases of 
adiabatic expansion. Assuming an equilibrium distribution fff misconceives the evolution 
of the temperature in the sense that higher temperatures and thus more high-momentum 
neutrinos are expected at earlier times. By contrast, the actual distribution takes the 
temperature plateau into account and is, hence, aware of the correspondingly less efficient 
production at high momenta. 
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# 


a 




ctr 




CR 


[M x ] 


1 


27 


8.4 x 10~ 4 


1.0 x 10" 


-3 


1.2 


7.0 


2 


210 


1.0 x 10~ 2 


1.4 x 10~ 


-2 


1.3 


5.3 


3 


1500 


0.29 


0.25 




0.87 


4.4 


4 


3500 


0.75 


0.24 




0.32 


3.4 


5 


6500 


1.0 


4.0 x 10~ 


-2 


3.9 x 10~ 2 


2.5 


6 


12000 


1.0 


4.0 x 10" 


-5 


3.9 x 10~ 5 


1.7 



Table 5: Indicators for the deviation of the thermal neutrinos from kinetic equilibrium at different times: 
M-R; a R aiJ d are introduced in Eqs. (B.3) and (B.5), cr is defined as cr — <jr//ir. The values of the 

T^and FigM 



scale factor correspond to a RH , a\, a%, 04, ag and an, cf. Section 



C Reheating temperature 

The energy transfer to the thermal bath, i.e. the reheating of the universe, becomes 
fully efficient when the nonthermally produced N\ neutrinos decay into standard model 
radiation. This happens once the Hubble rate H has dropped to the value of the effective 
decay rate rjy- of the nonthermal Ni neutrinos. The temperature at this time, t = tjui, 
defines the reheating temperature Trh 



RH, 



H(t 



RH) 



T, 



RH 



T(t 



RH) 



(C.l) 



Notice that this coincides with the common definition of Trh in scenarios in which the 
universe is reheated through the decay of some species with effective decay rate T. 

The decay rate relevant to our scenario, rjy- , corresponds to the T = neutrino 
decay width weighted with the average inverse time dilatation factor for nonthermal 



neutrinos (cf. Section 3.2.3) 



ATi 



It 



Em 



7 = lt(t RH ) ■ 



(C.2) 



In the course of our numerical analysis Eqs. (C.l) and (C.2) are used to determine the 
reheating temperature as a function of the model parameters, Trh = Trh(vb-l, Ml, mi). 
The result of this computation is presented in Fig. [5] in Section |5j In this appendix we 
shall illustrate how one can understand the exact outcome of the Boltzmann equations in 
terms of increasingly accurate approximations (also shown in Fig. [5]). 



We begin with Eq. (C.l ) and try to solve it analytically for Trh- First, the Friedmann 
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equation allows us to express Hit^u) through the total energy density at time tjui 

tRH 



RH 



3M p 2 



Ptot(t 



RH 



(C.3) 



Here /3 is a correction factor that accounts for the fact that we do not determine the scale 
factor a dynamically but simply approximate it by means of constant effective coefficients 
to in the equation of state. This imprecision in a is then transmitted to H such that 
it does not fulfill the Friedmann equation exactly. Second, let us denote the fraction of 
the total energy density that is stored in radiaton at time by a" 1 . With the aid of 



Eq. (27) we may then write 



H 2 (t 



RH) 



87T 7T 

3iW|30 



Combining Eqs. (C.l), (C.2) and (C.4) one finds 



RH 



a -l/4 /3 -l/2 7 -l/2 



90 



9*,p Trh ■ 



1/4 



(C.4) 



(C.5) 



By construction this formula yields the same results as Eq. (C.l ). Its main advantage over 



Eq. (C.l), however, is that it allows us to estimate Trh with varying precision depending 



on how accurately the correction factors a, and 7 are taken into account. 

When solving the Boltzmann equations numerically we have also determined these 
factors along the way. The dependence of the time dilatation factor 7 on the neutrino 



parameters is shown in Fig. 10 As all nonthermal neutrinos are produced with initial 



energy \ms — \f] 2 M% ~ 150Mi it is clear that 7 is bounded from above: 1 < 7 < 150. 
In practice, we find that 7 takes on values roughly between 1.1 and 88 entailing 7 -1 / 2 



factors in Eq. (C.5) approximately between 0.95 and 0.11. The general behaviour of 7 



as a function of the model parameters is mainly controlled by the ratio of the S and Ni 
decay widths 



7 



7(rV r *) 



r° 
1 5 



oc 



Mi 



(C.6) 



EW 



The larger compared to T° s the less contribute very long-lived nonthermal Ni neutrinos 
to 7. Most nonthermal neutrinos present at t = tun were then produced just shortly before 
and are hence relativistic. On the other hand, if the 5* bosons decay very fast, T° s 3> , 
the nonthermal neutrinos are mainly produced at the early stages of reheating and 7 is 
rather dominated by elder, nonrelativistic neutrinos. 
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Figure 10: Contour plots of the relativistic correction factor 7 as a function of the parameters rh\ and 
M\ for the three different choices of vb-l (cf. Eg. ( 102[ ) ). 7 _1 is defined as the average inverse time 
dilatation factor for nonthermal N\ neutrinos at t = t^u, cf. Eq. (C.2). The background colours reflect 
the parameter dependence of 7 as aiso indicated by the contour lines and labels. In the reddish regions 
(large 7 J the nonthermal neutrinos are relativistic at t = tun, in the bluish regions (small 7) rather 
nonrelativistic. In principal, 7 can take on values between 1 and |ms/-Mi ~ ^i] ~ 2 ~ 150. Its behavior 
in parameter space is mainly controlled by Tg/T^ oc (w|_ L m-i) / (y 2 EW Mi). 



55 



The correction factor a, the total-to-radiation energy density ratio at t = tun, in- 
creases when going to larger vb-l or fh% and decreases for smaller M\. It hence essentially 
shows the same trends in parameter space as 7. The physical reason for this is that large 
7 implies a rather long-lasting stage of N\ production through S decay which persists 
until shortly before the N\ neutrinos decay themselves. Thereby, the N\ abundance at 
t = ends up being still quite large. On top of that, in the case of highly relativistic 
neutrinos, the effective decay rate changes faster with time than for nonrelativistic 
neutrinos. Or, put into mathematical terms: jt (cf. Eq. (C.2)) is a monotically decreasing 
function of time which means that large values of 74 entail large values of ^7t _1 - Thus, 
if at t = tRH the decay rate Tff equals the Hubble rate H and it is fast changing it was 
much smaller before. Then not as many nonthermal neutrinos decay at times t < t^H and 
Pn^rh) may contribute much more to ptotiiRu) than Pr^rh)- Numerically, we find that 
a lies in the range between 3 and 4 for almost two thirds of the investigated parameter 
space. In approximately 85 % it is of O(10), in 10 % of (9(100) and in 5 % even larger 
up to a ~ 7.8 x 10 4 . In Eq. (C.5) the factor a -1 / 4 typically has a size between 0.75 and 
0.49. But in extreme cases it can become as small as a ~ 6.0 x 10 -2 . 

The correction factor in the Friedmann equation turns out to be quite constant in 
parameter space. We find that it varies between 0.53 and 0.92. Its standard deviation 
with respect to its mean value is rather small: /3 = 0.82 ± 0.06. A factor smaller than 
one is the expected consequence of our approach to the calculation of the scale factor: 
After t = ts we assume pure radiaton domination, u = p/p = 1/3, although for times 
t > ts surely still some nonrelativistic S bosons contribute to the total energy density. 
This leads to an overestimation of the speed at which the Hubble rate decreases and in 
Eq. (C.3) to a too small Hubble rate compared to the right-hand side of the equation. A 
correction factor of /3 -1 / 2 ~ 1.1 in Eq. (C.5) is, however, almost insignificant. 

In conclusion, we can say that 7 constitutes the largest correction, followed by a, the 
factor /3 can almost be neglected. This observation leads us to three increasingly accurate 
approximations for the reheating temperature T^, T R l jj and T R 2 ^, all of which are also 



shown in Fig. [5] First, we set a = /3 = 7 = 1. With g+ p = 915/4 and Eq. (12), we then 
have 



T ™ = ( ) \/ T k M P ( c - 7a ) 



90 



87r s g. 



1/4 



0.2 x /T%M p (C.7b) 

\ 1/2 / 
mi \ 1 Mi 



8xlo GeV lio^v [mm ■ (C7c) 
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This estimate can be improved by including the time dilatation factor 7, 



T« =7- 1/a lffi^0.2^ 1 M p . (C. 
Finally, we drop the assumption that at t — t R n the entire energy resides in radiation, 



T R % = «- 1/4 7" 1/2 T^ * 0.2 Jt%M p . (C.9) 



(2) 

The remaining difference between the outcome of the Boltzmann equations and T RH - is 
then quantified by (3 

T RH =(3- 1/2 T^, (CIO) 
log 10 = log 10 t|£ + A log 10 7W , A log 10 T RH ~ 0.04 . (C.ll) 

To conclude, let us apply the above formula? to the specific parameter example which 
we discussed in Section|4j For the parameter values listed in Tab. |2j solving the Boltzmann 



equations leads to a reheating temperature of T RH = 4.1 x 10 GeV (cf. Eq. (97)). This 
result can be compared to the three estimates introduced in this appendix. For the three 
correction factors a, and 7 we obtain 

a -3.2, /3-0.84, 7 -14. (C.12) 

At t — t R H there is, hence, still roughly twice as much energy in nonthermal neutrinos as 
in radiation. Equality in terms of the energy content is not reached before a ~ 4000 in 
our parameter example. Moreover, with a typical energy ~ O(10) Mi the nonther- 



mal neutrinos are clearly still relativistic during reheating. According to Eq. (C.12) the 
estimates T R ^, and T R 2 ^ yield reheating temperatures of 

T R % ~ 1.9 x 10 10 GeV , T R 1] H ~ 5.0 x 10 9 GeV , T R % ~ 3.7 x 10 9 GeV . (C.13) 

Notice that T R l multiplied by /3~ 1//2 reproduces again the numerical result for the reheat- 



ing temperature T R h in Eq. (97). We conclude that the most naive estimate T RH - is off 
the actual value by roughly an order of magnitude. T R jj and T RH - respectively deviate 
from Trh by 10% and 20%. 

D Semi-analytic reconstructions 

Our study of the parameter space in Section [5] allowed us to determine the Ni neutrino 
mass Mi and the reheating temperature T R h as functions of rhi and mg such that the 
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gravitino abundance always has the right size to account for dark matter (cf. Figs. [7] and 
8|. In this appendix we now attempt to reconstruct these results by means of simple 
analytic expressions and with the aid of our numerical findings for tjb and T R h- 

As gravitinos are nonrelativistic, their present contribution to the energy density of 
the universe is given as 



Q d h 2 



tt d h 2 (v B _ L ,M 1 ,m 1 ,m d ,m~ g ) = mg^gn° h 2 /p c 



n|/n°. (D.l) 



In order to relate the gravitino-to-photon ratio r] G to the corresponding number densities 
during reheating we make two simplifying assumptions. First, we say that after t = t R u 
the entropy of the thermal bath is not increased much further which leads us to 



6x 



a 



-Q-ny{tRH) ■ 



(D.2) 



Second, we assume that at t — t R n the gravitino production becomes inefficient such that 
at later times not many further gravitinos are produced, 

3 



4 



8, 



a (t 



RH, 



«0 



n G (t 



(D.3) 



Meanwhile, this second assumption also implies that at t — t R u the gravitino production 
rate 7g is of the same order as the Hubble rate H 



7g(< 



RH) 



8 3 l H{t RH ) & n 5 (t RH ) = 8- 



7g(* 



RH) 



H{t 



RH j 



(DA) 



n G (t RH ) 

The three correction factors 8i > 1, 82 > 1 and 83 ~ 0(1), introduced in Eqs. (D.2), 
(D.3) and (D.4), respectively, quantify the deviations of the actual values of n°, n°~ and 
n c(^RH) from our approximations. Combining them in one factor 8 = 6283/81 we may 
write for r\ G 

" 7g(*i 



8 



<^RH) 



(D.5) 



where n 7 (taff), l G [tRH) and H{t R u) directly follow from Eqs. (27), (86) and (C.4). In- 
serting Eq. (D.5) back into Eq. ( [D.l ) we find for VL Q h 2 



n d h 2 



zf G ( T RH) 



m G + 



ml(T RH ) 
2>m G 



T, 



RH 



a 



- i / 2 r 1 5. 



(D.6) 



where f G (T RH ) stands for 



f G ( T RH) 



90 



1/2 



^g 2 s {T RH ) 
gjM p 



In 



T 2 

1 RH 



m 2 (Tj 



RH, 



+ 0.8846 



(D.7) 
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Eq. (D.6) may conveniently be rewritten as 

Trh 



Sl 5 h 2 = eC^Trh) 



10 9 GeV 



C 2 (Tj 



RH 



m 5 



10 GeV 



10GeV\ / m g (u, ) 
^TJ V800 GeV 



(D.8) 



with C\ and C 2 being denned as 



C2\Trh) 



(800 GeV) 1 
10 GeV 
10 GeV 

(800 GeV) 2 



x 10 9 GeV x 



dpi 



RH 



x 10 GeV x 



x /fi(T, 



(D.9) 
(D.10) 



The dependence of G\ and C% on the reheating temperature is presented in Fig. 11 



We find that Ci/C^ ~ 0(100) which means that for nig ^> mg the term linear in in 



Eq. (D.8) can usually be neglected. Notice that doing so and setting e — 1 turns Eq. (D.8) 



into Eq. ^ in the introduction. 



Confronting Eq. (D.8) with our numerical data shows that e usually differs from 1 and 



is slightly parameter-dependent, preventing us from determining it a priori. However, e 
can be determined a posteriori. In the region of parameter space in which fl^h 2 = QpMh 2 
we find that e = 1.20 ± 0.24. Restricting ourselves further to parameter values for which 
gravitino dark matter also is in accordance with successful leptogenesis we obtain e = 
1.32 ± 0.15. In Eq. ^ such a correction factor would be reflected in a change of the 
numerical coefficient from 0.26 to 0.34. 

Eq. ( D.8[ ) then implicitly determines the reheating temperature as a function of 



and ra<j(/xo)- Fixing the gluino mass at 800 GeV and solving Eq. (D.8) for Trh provides 



us with an estimate for the reheating temperature solely dependent on mg (cf. Fig. 12). 
This is all we need to be able to reconstruct Figs. [5] and [6j With Trh = TR H (m^) at 



hand we can compute the reheating temperature for all values of the parameter triple 
[y b—l ; f^i ; f^Q) ■ From our numerical results for Trh as a function of vr-l, Mi and fhi, 
shown in Fig. [5j we then infer the corresponding values of M\. Our results for i]b in Fig. [6] 
finally point us to the respective baryon asymmetries 



Q d h 2 



0.11 



Trh = Tjuj(jtIq) , 

Mi — Mi (vb-l: Trh, fhi) 

Vb = Vb (v B -L,Mi,mi) . 



(Eq. flPJD ) 

(Fig.§ 

(Fig.g 



The outcome of this procedure is presented in Figs. 13 and 14 As it turns out our 
reconstructed results fit the exact data from the Boltzmann equations amazingly well. 
We thus conclude that our numerical results for the reheating temperature Trh and the 
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Figure 11: Dependence of the coefficients C\ and C2 in Eq. (D.8) on the reheating temperature Trh 
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Figure 12: Analytic estimate for the reheating temperature Trh as implicitly dchncd by Eq. (D.8) with 
e = 1.32 for varying gravitino mass mg and fixed gluino mass, mg = 800 GeV. 



60 



baryon asymmetry 7]b when combined with Eq. (D.8) essentially suffice to reproduce 
Figs. [7] and [8} 



Let us check how well the result for Vt^h (cf. Eq. (D.8)) that we obtained for the 



parameter example discussed in Section [4] can be reproduced with the formulae developed 
in this appendix. The three correction factors Si, 82 and $3 turn out to be 



Si 



N 7 (t 



2.7, So 



RH, 



15 , 5, 



H(t RH )nn{t RH ) 



7g(* 



RH , 



0.35, (D.ll) 



which implies that entropy production has almost completed at t — t R H- By contrast the 
gravitino production rate is still roughly three times as large as the Hubble rate at this 
time such that the bulk part of the gravitinos is, in fact, produced at the last stages of 



reheating and later. In combination with Eq. (C.12) the three factors yield values of 5 
and e of 



S = S 2 S 3 /S 1 ~ 2.0 , e = a- 1/2 (3- l 5 ~ 1.3. 



(D.12) 



This result for e coincides with the fit value used for the reconstruction of Figs. 13 and 
14 Notice also that the effects of the various approximations parametrized by a and S 
tend to cancel such that overall factor e represents a correction of only 30% in the end. 



Given the reheating temperature in Eq. (97) the coefficients C\ and C 2 in Eq. (D.8) take 
on the following values 



Ci{Z 



RH) 



0.17. 



C 2 {T RH ) ~ 2.0 x 10" 



(D.13) 



Based on Eq. (|D.8|) we can then estimate the gravitino abundance 

l()GeV\ / m^p) V 



Q d h z = 1.3 x 0.17 



T, 



RH 



10 9 GeV 



0.002 



m G 



10 GeV 



+ 



m 5 J \ 800 GeV J 



1.3 x 0.17 x 4.1 x (0.02 + 0.1) ~ 0.11 



(D.14a) 
(D.14b) 



which is exactly the value we obtained solving the Boltzmann equations. Without the 
correction factor e, we would end up with a too small value, flgh 2 ~ 8.5 x 10~ 2 . If we were 
to neglect the term linear in mg in addition, our estimate would come out even smaller, 
7.1 x 10~ 2 . In concluding, we also mention that our result for the reheating 



tt 5 h 2 



temperature in Eq. (97) coincides, by construction, with the value required for gravitino 

100 GeV corresponds a temperature of 

Trh ~ 



dark matter: In Fig. 12 we read off that to mg 
4.0 x 10 9 GeV. 
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Figure 13: Semi-analytical reconstruction of Fig.^(solid lines) on the basisofEq. (D.8) withe = 1.32 and 
the numerical results for Trh and t]b- For comparison also the Mi contours deduced from the Boltzmann 
equations (dashed lines) are shown. They deviate from the reconstructed results as the correction factor 
e = a~ 1 / 2 /3~ 1 S actually is parameter-dependent. 
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Figure 14: Semi-analytical reconstruction of Fig. [$] (solid lines) on the basis of Eq. (D.8) with e — 1.32 
and the numerical results for Trh and t/b- For comparison also the Trh contours deduced from the 
Boltzmann equations (dashed lines) are shown. They deviate from the reconstructed results as the 
correction factor e = or 1 ! 2 'f3~ 1 5 actually is parameter-dependent. 
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Finally, our results may be easily generalized to gluino masses other than 800 GeV. In 
fact, for given values of vb-l, fn>i, Mi and it is possible to keep r\ B and fl^h 2 constant, 
when changing rrig, by simply rescaling the gravitino mass, 



m~ — > 



m,Q (m,,. in 



"gi ""a 



s) 



m d (800 GeV, m~) 



(D.15) 



As for the baryon asymmetry, this is a trivial consequence of the fact that 7]b is a function 
of Vb-l, fn>i and Mi only. In the case of the gravitino abundance we observe that for 
fixed reheating temperature, Trh = Trh (vb-l, fn>i, Mi), Q^h 2 remains constant as long 
as mg(mj,m~) is chosen such that the term in square brackets in Eq. (D.8) does not 
change, 
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m d J V 800 GeV 



(D.16) 



where Ci = C*2(Tr#) is defined in Eq. (D.10). From Eq. (D.16) we can determine the 



rescaled gravitino mass as a function of the rescaled gluino mass nig and the original 



gravitino mass m~. As Eq. (D.16) is a quadratic equation in mg, it generically has two 



solutions m~, one of which is typically closer to the original gravitino mass than the 
other. m°~ lies right in between m~~ and mX once the two terms in square brackets in 

G '■" G G 



Eq. (D.8) are of equal size, i.e. when gravitinos in helicity ±| states contribute exactly 



as much to the total abundance as gravitinos in helicity ±| states. One easily sees that 
this is the case when m~ ~ 230 GeV (cf. Eq. (117)). When going to larger trig, we have 



m % ~ m % ^ m p=, above 230 GeV and m% < m~ m~ below 230 GeV. At m„ smaller 

Cr G G G G G 

than 800 GeV, we always find m~ < m°~ < m~. 

G G G 

If the gravitino mass is much smaller than the gluino mass, almost only the goldstino 



part of the gravitino is produced and the term linear in mg in Eq. (D.8) can be neglected. 
The scaling behaviour of the gravitino mass then becomes trivial 



uiq <C m- g 



o / m g 



(D.17) 



Actually, the rescaled gravitino mass mg also is a function of T RH as it depends on 
the coefficient Cz^Trh)- But as discussed in this appendix, there is an almost unique 
correspondence between the gravitino mass and the reheating temperature. In order to 



solve Eq. (D.16) we thus simply read off Trh from Fig. 12 as a function of the input 
gravitino mass, T^h — Trh ( m %) ■ As can be seen from Fig. 14 , this simplified reheating 



temperature is in good agreement with the exact outcome of the Boltzmann equations. 



Our solutions m~ for the rescaled gravitino mass are presented in the two panels of Fig. 

G 



15 
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Figure 15: Contour plots of the two solutions (mi-) of Eq. (|D.16|) for the rescaled gravitino mass m- ( 



a function of the rescaled gluino mass wig and the original gravitino mass m~. The black solid contours 
correspond to constant values of mg (given next to the green dots). They serve as level curves that 
allow a determination of for arbitrary points in the (m^m-j -plane. They can also be regarded 
as function graphs of mP~ as a function of nig for constant m^. We restrict ourselves to the interval 



20 GeV < m ~ < 700 GeV in this figure. Below 20 GeV, Eq. ( D. 17 1 provides an excellent approximation. 



respectively. In the gray shaded regions there are either no real solutions of Eq. (D.16) or 
the rescaled gravitino mass is larger than the corresponding gluino mass, > m- g . The 
former case implies that it is impossible to keep the gravitino abundance constant when 
going to larger m-g while sticking to the reheating temperature TR#(m~). In the latter 
case, the gravitino would not be the LSP any longer. 
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